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The paper' investigates optimum processes in systems whose 
behavior is described by difference boundary problems for partial 
differential equations. 


The majority of physical processes which are encountered by the 
engineers during their practical activity can be controlled; and, 
consequently, during their realization, one strives to attain the optimum 
(in some sense) alternative. The maximum principle of L. S. 
Pontryagin 1 appears as a mathematical method for the solution of 
problems of optimum control when the processes may be described by 
the ordinary differential equations. However, numerous control 
processes are described by partial differential equations with additional 
(boundary of initial) conditions. These equations may be of diverse 
type (equations of mass or heat exchange, equations of hydro or aero- 
dynamics, heat transfer, kinetics of chemical reactions, etc). If the 
behavior of the control system is described by equations among which 
there are some with partial derivative^, then it is called a system 
with distributed parameters 2 . In numerous simpler cases, such 
systems may be described by differ ential -difference s equation; and, 
consequently, one can still apply the maximum principle . 


The problems of optimum control of more complicated systems 
can not be solved directly by means of the maximum principle of L. S. 
Pontryagin (see article by Kharatishvili 4 , pp. 5 16-5)8). Consequently, 
attempts have been made to generalize this principle in such a way that 
by its application one could investigate control processes of more 
complicated systems with distributed parameters 5,6,7 ’ 8 ' 9,10 ’ 11,12,13,14,1 ". 
In particular, the work 10 proposed a method based on the use of 
differential equations in Banach spaces. In numerous cases, such an 
approach allows the investigation of partial differential equations as 
if they were ordinary differential equations and solves the problem of 
optimum control using as the optimality criterion the functional 


I = J f(t, X (t), u(t)) dt 
p 


( 1 ) 


In spite of obvious advantages, this method has also substantial short- 
comings since the introduction of Banach spaces requires additional 
limitations on the class of permissible controls which are not caused 
by the physical essence of the problem. In addition, the choice of the 


'‘'The basic content of the work was presented at the seminar of 
L. S. Pontryagin on theory of optimum processes , 13 February 1964. 
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functional (.1) as the optimality criterion for the problem with partial 
differential equations is not as successful as in the case of problems 
with ordinary differential equations. In particular, the indicated 
method does not solve the important practical problem in which the 
functional ( 1) is substituted by an integral evaluated over a surface 
bounding the domain with which the equations are under investigation, 

1 3 

Of definite interest is the method based on the representation of 
control quantities by means of integral relationships. However, its 
justification can not be viewed as satisfactory. In addition, the applica- 
tion of this method to processes described by boundary problems with 
partial differential equations can not be considered sufficiently efficient 
because of the following reasons. First, the reduction of boundary 
problems to integral equations can not always be carried out though the 
problem can still be solvable by other methods. Second, it is always 
desirable to have optimality conditions expressed directly through 
quantities entering the equations and additional conditions. 

In the present work, we used the method of solution which is to 
an equal degree applicable to the cases of hyperbolic, parabolic, and 
elliptic equations. Using this method, L. I. Rozopoer 16 studied the 
case when the control process is described by ordinary differential 
equation and finite differences equations. In subsequent works 17 , we 
obtained invariance coefficients for systems relative to the external 
interactions whereas the starting point for investigation was a formula 
for functional increments found in the work of Rozopoer 1 . Analogous 
results (though only for special cases) were obtained also for systems 
with distributed parameters. 

The paper consists of five sections. In Sections I and II we 
investigate various problems of optimum control of processes 
describable by boundary conditions for hyperbolic equations with data 
on characteristics. The necessary optimality conditions are formulated 
in the form of a maximum principle. 

In Section III we establish connection of the investigated problem 
with the problems of variational calculus. It is shown that the Euler- 
O strogradskiy equations can be derived from the maximum principle 
if the control domain coincides with the entire space. If this domain 
is closed, then along the optimum surface one may find that the 
Legendre conditions are even not satisfied. 
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In Section IV we study problems of optimum control when processes 
are described by the boundary conditions for parabolic systems. We 
obtain a formula for the functional increment by means of which one 
finds the optimality conditions. These results are extended to analogous 
problems connected with the elliptical and hyperbolic systems. 

Section V is devoted to the problems of the invariance theory. For 
linear equations we obtain the necessary and sufficient conditions for 
invariance relative to the external interaction whereas the criteria of 
invariance we choose functionals analogous to those investigated in 
Sections I through IV. 

The author uses the occasion to express his gratitude to L. S. 
Pontryagin and the participants of his seminar for their interest in the 
present work. In addition, the author is sincerely grateful to V. G. 
Boltyanskiy, O. A. Oleynik, and Yu. V. Yegorov for very useful 
discussions of the results obtained in the paper. 
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Section I. OPTIMUM PROCESSES IN SYSTEMS WHOSE BEHAVIOR 
IS DESCRIBED BY HYPERBOLIC EQUATIONS 


1. Formulation of the Problem. Optimality Conditions 

Let the control process be described by a system of equations 

z ixy “ ^i ( x > y> z i » » z m» z ix> ••• » z mx> Z iy’ * ** ’ z my» v )> 

i = 1 , ... , m, (1.1) 

where the functions f^ have within the domain G (0 < x < X, 0 < y < Y) 
continuous derivatives of the first order in x and y and twice continuously 
differentiable over other arguments. As a class of permissible controls, 
we use the set of sectionally continuous functions v v(x, y) defined 
within the region G and with values within a certain bounded convex 
domain V (open or closed) of the r - di men sional euclidian space. We 
assume that the line of di s continuity of the permissible control is 
sufficiently smooth. We impose onto the function defined by equations 
(1.1) certain boundary conditions (Goursat conditions) 

z i ( 0 , y) = m i(y) , z i(*, 0) « ^(x) - i - 1. ••• . rn, (1.2) 

where and are continuous, sectionally continuously differentiable 
functions defined within the domain G and satisfying matching conditions 

^i(O) = 4^1 (0) . 

Each permissible control may be associated with a unique solution 
z (x, y) = jzj (x, y) z m( x > V)\ 

of the problem (]. l)-(l.Z) having derivatives z* integrable over the 

is ^ y 

domain G (see work of Budak and Gorbunov ). However, here one 
should distinguish two cases. 

1) If the line of discontinuity of the function v(x, y) is parallel 
to one o f the coor di nate axes , the boundary problem ( 1 * l)-(l.Z) 
decomposes into two analogous problems within the regions bordering 
one to another along that line. By solving consecutively these problems, 
we determine the solution of the original problem which may be 
continuous within the domain G and everywhere except along the points 
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of the discontinuity line of control v(x, y) , and will have continuous 
derivatives z^ x (x, y), z^y (x, y), and zq X y(x, y) (see work of Budak 
and Gorbunov 1 ). 

Z) Let the discontinuity line r of the function v(x, y) have no 
common sections with the characteristics of the system (1. 1) over 
any segment different from zero. Under the solution of the boundary 
problem (1. l)-(l.Z) we understand the function z (x, y) which satisfies 
the system of equations (1. 1) in all points of the domain G not lying on 
r and the conditions (1.2) together with certain a priori given conditions 
of smoothness along r (see work of Yegorov ). Such a solution is 
uniquely determined; it is continuous within the domain G and has 
sectionally continuous derivatives Zj x , z^y, and Zj X y. 

Consequently , in what follows we will assume that to each per- 
missible control one can associate a class of functions for which the 
boundary condition (1. l)-(l.Z) can be uniquely solved. 

Let A], i = 1, ... , m, be a given system of real numbers. Let 
us take an arbitrary permissible control v(x, y) and denote by z (x, y) 
the respective solution of the problem (1. l)-(l.Z) and study the 
functi onal 


m 

S =£ A A /-i (X , Y) , 

i = 1 


( 1 - 3 ) 


where X and Y are constants entering into the definition of the domain G. 

The problem is then: among all the permissible controls one is 

supposed to find such a control v(x, y) (if it exists) that over the solution 
z(x, y) of the Goursat problem corresponding to this control the 
functional S attains its largest (smallest) value. 


We will call the permissible control realizing the minimum 
(maximum) of S the min-optimal (max- optimal) control according to 
S (see work of Rozopoor 16 ). 


Note that the problem ( 1 . 1 ) - ( I . Z) under consideration is of great 
theoretical and practical interest. The study of the solvability of this 
problem for various assumptions relative to the function f^, ^ , and 4 d 
is the object of extensive 1 iterate r e 20 ’ 21, 22> 23, 24 ’ . It is also known 26,27,28 

that the study of sorption and desorption of gases, drying processes, 
and the like reduces to such a problem. The presence of control para- 
meters in Equations (1.1) allows the control process, and in numerous 
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cases the selection of the best operating conditions which (from the 
mathematical point of view) reduces to the calculation of the maximum 
or minimum of a certain function. In numerous cases the problem may 
be reduced to the study of the functional (1. 3). Let us investigate certain 
example s . 

1) One is supposed to minimize the functional 

1 = // f o ( x > y> z > z x> z y> v ) dxdy . 

G 

If we introduce a new variable z Q by putting 

z oxy = y > z > Z X. Zy. V), z 0 (0, y) = z 0 (x, 0) = 0 , (1.4) 


the problem reduces to the calculation of the minimum of the function 
S = z Q (X , Y) which represents a special case of the functional (1.3) 
and is defined over the functions z Q , ... , z m specified by the totality 
of the relationships (1. l)-(l.Z) and (1.4). 


Z) One is supposed to minimize the functional I = ^(z 1 (X, Y), 

• * • » z m (X > Y)), where $ is a twice continuously differentiable function. 

We introduce a new function z Q (x, y) by means of the equation 

y), ... , Z m (x, y)) 


m 


J Oxy 


E 

i , k = l 


9 z i9zk 


z ix z ky 


m 

V' 9$ 

+ Z-/ fi(x ’ y> z> z x> z y’ v) 

i = 1 


and additional conditions 


z 0 (°> y) = (y) ®m(y))> z 0 ( x > °) = $ (4- , 1 ( x ). .... ^ m ( x ))- 


This reduces the problem to the study of the functional S = z 0 (X, Y). 



3) One is supposed to minimize the functional 
X 

I = f F (x, z(x, Y), z x (x, Y)) dx . 
o 

We introduce the auxiliary function z Q (x, y) by means of the 
equation 


! 


z oyz 




8F 

8 Z-: 


£i(x, y, 


IX 



and additional conditions 


x 

z 0 (0, y) - z 0 ( x > °) = / F ( x > + ( x )> 4*' ( x )) dx . 

0 

4) In an analogous manner one investigates the problem 
concerning the minimization of the functional 

Y 

I = / F (y , z (X , y), z (X, y)) dy . 

o 

For the solution of the formulated optimum problem one introduces 
the auxiliary functions Uj , ... , u m by means of the equation 


u ixy 


9H(x, y, p, v) 
9zj 


d / 9H(x, y, p, v) 
dx \ 9z ix 


d / 9H ( x , y, p, v) 
dy\ 3 z iy 


(1. 5) 


and additional conditions 


u ix ( x ’ Y ) _ - 


9H (x, Y, p, v) 


9z- 


L y 


^i y ( X > y) 


9H (X , y, p, r) 

9 z ir 


, u i (X , Y) = Ai 


( 1 . 6 ) 
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where are constants entering into the definition of the functional S, 

P " ( z 1 > • • • ) z m> a i > • • • * u m ’ z i x ’ * * * 9 z mx » z j y > * • * » z my ) * 

H = 2 u i f i ( x . y. z - z x > z y> v ) 

Conditions (1.6) represent ordinary linear differential equations 
with initial conditions. 

In the general case, the right-hand side of the equations (1. 5) 
contains functions zq xx , z^yy, v x , and Vy. However, from the conditions 
imposed on Equations (1. 1) and the permissible controls, it does not 
follow that such derivatives must exist. Consequently , we will assume 
in what follows that the function f^ may be represented in the form 

m m 

f i 2 a ijk( x > y> z > z kx z jy + X) b ij (x ’ y> z) z j x 

j, k = l j = l 

m 

+ ^ c i j(x, y, z) Zj y + (x, y, z, v) . 

j " i 

where the functions a — b — , c - , and d^ are continuously differentiable 
over x and y and are twice continuously differentiable over the other 
arguments. If it turns out that a — by , and c— do depend on v, then 
one must require that the permissible controls have sectionally 
continuous derivatives v x (x, y) and Vy(x, y). 

While satisfying these conditions, the system of linear equations 
(1.5) has sectionally continuous coefficients and together with the 
additional conditions (1.6) determines uniquely u x (x, y) , ... , u m (x, y) 
for each permissible control. Consequently, in what follows we will 
assume that the function f^ and the permissible controls are such that 
the boundary problem ( 1 . 5) - ( 1 . 6) be unique solvable for each per mis sible 
control. 

We will say that the permissible control v(x, y) satisfies the 
condition of maximum if the relationships 

H (x, y, p(x, y), v(x, y)) ((=)) sup H (x, y, p(x, y), v) , (1.7) 

v €IV 
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is satisfied with z(x, y) and u(x, y) - the solutions of problems (1. 1)- 
(1.2) and ( 1 . 5) -( 1 . 6) cor r e s ponding to the control v (x , y ) , whil e the 
symbol ((-)) denotes equality valid in all points of the domain G 
(0 < x < X, 0 5 y < Y) while there may exist sets of points lying over 
a finite number of lines with zero surface. The conditions of minimum 
is defined in an analogous way. 

THEOREM 1 (the principle of maximum). For the permissible 
control b(x, y) to be min-optimal (max - optimal) according to S, it is 
necessary that it satisfies the condition of maximum (minimum). 

Although this theorem does not supply sufficient conditions for 
the existence of optimum control, it maybe utilized for the practical 
solution of the optimum problem. As a matter of fact, the solution 
of this problem, according to the principle of maximum, leads to the 
need of determining 2n 4 1 unknowns z{ , uj_. and v from the 2n 4 1 

Equations (1.1), (1.5), and (1.7). First 2n relationships represent 
second order differential equations whose solution in general generates 
4n arbitrary functions. To eliminate them we have 4n additional 
conditions (1.2) and (1.6). This is sufficient to define, generally 
speaking, the separation of the solution of the problem (1. ])-( 1. 2) 
satisfying the conditions of the maximum principle. If it appears from 
the meaning of the problem .that the optimum problem must have a 
mandatory solution, then at least one of the discovered isolated solutions 
must be the desired one. 


2. The Formula for the Increments of the Functional S 

To prove Theorem 1, we study the functional 

Mp. 'i - // 

G 

If v is a certain permissible control and z = z (x, y) a solution of 
the problem (1. l)-( J.2) corresponding to this control, then the functional 
I is equal to zero for an arbitrary function u = (vij , ... , u m ). 

Let v = v (x, y) he a certain permissible control and z(x, y) and 
u (x , y) be the solution of the boundary conditions ( 1 . 1 ) - ( 1 . 2) and ( 1 . 5) - 
(1,6) corresponding to this control. Let us give the function v an 
arbitrary permissible increment Av and let us denote by z 4 Az and 
u 4 Au the solution of the same problems but corresponding to the 
control v 4 Av. It is clear that the functions A z[ and Aiq satisfy the 
conditions 


u; z; 


ixy 


II (x . y, p, v) 


dx dy 
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A A 8H 

A z- = A 

1X Y DU: 


A u jxy " A 


9H _d_ / 9H \ _ _d_ / 

9zi dx \ 9z ix / dy \ 


and the additional conditions 


( 1 . 8 ) 


Azj (0, y) = A (x , 0) = 0 


(1.9) 


Au ix (x, Y) = 


A 9H(X ’ z Y ’ P ’ V) . Au (X, y) 

9 z i y 


A Uj (X , Y) = 0, i = 1 


9H (X, y, p, 


(1. 10) 


(1. 11) 


where 


8H 9H (x, y, p + Ap, v + Av) 9H(x, y, p, v) 


( 1 - 12 ) 


Equations (1. 10) are ordinary differential equations where with 


linear f- functions 


Au-(x, Y) = 0, A Uj (X , y) = 0 


( 1 . 13 ) 


are their solutions which satisfy the additional conditions (1. 11). 
Because of the uniqueness theorem, functions (1. 13) form a unique 
solution of the boundary problem (1. 10) - ( 1. 11). 

Consequently, in accordance with the above remark, 


AI = I[p + Ap, v + Av] - I [ p, v] = 0 


( 1 . 14 ) 


On the other hand, 


A 1 ^ | A uq A + uj. A Zq X y 3 Aq z ixy| 


H ( x , y , p + A p , v + A v ) - H ( x , y , p , 


„ I l 


v; I s 


dx dy . ( 1 . 15) 
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The expression under the sign of the integral is transformed by 
means of the Green's formula (see work ofTrikomi 29 , p. 196): 

// ( < 9 s xy - s q X y) dxd Y = / (q s y - s 9y) d y - (qs x - sq x ) dx , 
G L 


where L is the contour limiting the domain G, and q and s - arbitrary 
functions having sectionally continuous derivatives of the first and 
second order. Since G is a rectangle, the Green's formula may be 
reduced to the form 



X Y 

/ Y r x 

(sq x )y = 0 dx -J (sq )* = 0 d y • ( 1 . 16 ) 

o v ' o ' y 


Let us insert into this equality q = Au^ and s = Azj, Taking into 
account Equations (1.8) with additional conditions (1.9), (1.10), and 
(1. 11), we obtain after elementary transformations 


m m 

J J X Au i A z ixy ^ d Y = // E 
G 1 = 1 G i = i 



A z£ 


+ A 


8H 

8zi x 


A z 


IX 


+ A 


3H 




ly 


J iy 


dx dy . 


On the other hand, because of the first m equations of (1.8), we have: 

. . rci . . rn ot t 

// E Au i Az ixy dxdy =// £ A — A ui dx dy . 

G i = 1 G i = 1 1 

From the two last equations we obtain 

4m 

// E Au i Az ixy dxdy = y // E A A Pi dxdy . (1. 17) 

G i - i G 1 = 1 Pl 

We now substitute into Equation (1. 16) q = u^ and s - Azj_. Then 
because of Equations (1.5) and ( 1 . 8) and the boundary condi tions ( 1 . 6) 
and (1.9) we have 



Since the function forms the solution of the system of equations 
( 1 . 1 ) , we have 

™ rr m 

JJ L, Au iZi dxdy = JJ Y, T— Au i dxdy . (l.K 

G i = l G i = i 1 * * 

Using the Taylor formula, we obtain the equality 
H(x, y, p + Ap, v + Av) - H(x, y, p, v) 


8H (x, y, p , v + Av) 

2 8^ Ap i 

1—1 11 


i 4m 

+ T E 


9 H (x, y, p 4 0A p, v 4 A v) 


i, k = l 


3Pi 9 Pk 


A pi A p k 


4 H (x, y, p, v 4 A v) - H(x, y, p, v), 0 < 9 < 1 . 


(1. 2C 


From Equations ( 1 . 14) , (1.15), and ( 1 . 1 7)-( 1 . 20) , it follows that 


m 


AI = - E A i Az i (X, Y) -ff [H (x , y, p, v 4 Av) 

G 


1 = 1 


i - 4 m n 

- H (x , y, p, v)] dxdy + j Jj E ' 

G i = i 


8H (x, y, p + Ap, v + Av) 


3Pi 


9H(x, y, p, v 4 Av)] 


f 9H ( x , y , p , v + Av) 



9H(x, y, p, v) 
api 


4 m 


A Pi dxdy - - // 2 


G'i, k = i 


9 H (x, y, p + 9A p, v + Av) 

a Pi a Pk 


Ap i Ap k dxdy * 


Applying to trie functionals 9H/9Pi the Taylor formula and taking 
into account the equality (1. 14), we finally obtain 

AS = // t H < x - Y> P. v + Av) - H(x, y, p, v)] dx dy - r\ (1.21) 

G 


where 

AS = £ A i A Zi (X, Y) 

is the increment of the functional S , h = q + r\ 

1 2 

A Pi dx dy , 


fs// 

i-i r. 


8H(x, y, p, v + Av) 9H (x, y, p, v) 


9pi 


3 Pi 


4m 


e /; 

i,k = l G 


9 H(x, y, p + 9Ap, v + Av) 9 2 H(x, y, p + 9j Ap, v + Av) 


3Pi3p k 


9Pi 3Pk 


Ap i A Pk dx d Y 


( 1 . 22 ) 


3. The Estimate of the Residual Term r\ in the Formula (1.21) 

To establish the necessary estimates of the quantity r\ we 
introduce auxiliary functions a^(x, y) and (3^(x, y), putting 

a i "" ^ z ix ’ Pi ~ ^ z iy 

Since the function T satisfies the Lipschitz condition, then from the 
first m equations of the system (1.8) and the condition (1.9) we obtain 


13 



I a i I £ N J E (|A Zi | + | Qi | + |0i|)dy + N, J ^ |Av k | dy , 

0 1 = 1 0 k = 1 

X 

|A Zi |< / | ai | dx , 

0 

(1.23) 

x m x r 

I ®i| " n/Z <l A *i| + KI + I 0 il> ** + N J X |Av k |dx , 

o i — i o k = l 

y 

l Az il- / l B i! d y > 

0 

where N and Ni are definite positive constants. Introducing the notation 
m m m 

a = E I Q il ’ P = E l 8 il ’ Y = E l Az il • 

i = l i = i i = l 

r 

Av = E l Av kl > (1.24) 

k = l 


we obtain from inequalities ( 1. 23) 

y p 

a (x , y) < NmJ* a(x, y) dy + NmJ [P(x, r|) + y(x, r|)] dr| 

o o 


+ N j m J* A vdy , y < f a (x, y) dx 

0 o 

(1.25) 

x § 

(3 (x, y) < Nm f (3 (x, y) dx + NmJ* [a (x, y) + (3(x, y)] dx 
0 0 

x X 

+ Nj mj* A vdx , y < f (3 (x, y) dy , 



where 


0 < x £ £ £ X, 0 £ y £ h £ Y . 


From this, because of the known lemma (see article by Nemytskiy an,d 

3 0 

Stepanov p. 19), it follows that 

h 


Y 


a(x, y) < M / [ y (x, y) + ( 3 (x, y)] dy + M ,/ Av(x, y) dy , 

o o 

| X 

(3 (x , y) £ P J [ Y ( x , y) + a (x, y)] dx + pJ A v (x, y) dx , 


where M, Mj, P, Pj are positive constants. Taking into account the 
estimates (1.25) for the function y> we obtain 


a (x , y) < M 2 J* (3(x, y) dy + mJ Av(x, y) dy , 

0 o 

% * 

P(x, y) < P 2 J a (x, y) dx + PjJ Av(x, y) dx . 


From this we find that 


il f Y X 

q(£, T|) < M 3 J f a(x, y)dx dy + M 4 j* A v ( x . y) dx dy 


o o 


Y 


0 0 


+ Ml J A V (£ , y) dy 
o 

t, e YX 

s(^,'n) - p >// 6 (x , y)dx dy + P 4 // A v(x, y) dx dy 


0 0 


0 0 


( 1 . 26 ) 


X 

+ Pi J Av (x, r \ ) dx 
0 

Integrating the first of these inequalities over £ between the limits of 0 
and £ and applying the above mentioned lemma, we obtain 
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J a(x, y) dx — M U J Av (x, y) dy dx 
o oo 

From this and the first inequality (1. 26) we have 

X Y Y 

a(x, y) - M 6 J* J* Av (x, y) dy dx + M 7 J* A v (x, y) dy, 
00 0 

0 < x < X, 0 — y — Y . 

In an analogous manner we find: 

X Y X 

3(x, y) < P 6 f f Av (x, y) dx dy + P 7 J Av (x, y) dx 
00 0 

From this and the inequality ( 1 . 25) we obtain 

X Y 

Y( X , y)<-Q f f A v (x , y) dy dx . 

0 0 

In this manner, because of (1.24), the inequalities 


|A zi (x, y)| - Q f f A v (x, y) dx dy , 

G 

Y 

|A z ix ( x > y)l - Qi If Av (x, y) dx dy + RlJ* Av (x, y) dy , ( 

G o 

X 

|Az iy (x, y)| < Q 2 // Av (x, y) dx dy + R zf Av (x 7 y) dx . 

G o 

are valid for all x and y(0^x~X, O^y^ Y). 

Applying an analogous approach to the last m equations of the 
system ( 1 . 8 ), we obtain 


Auf (x, y ) | ■o,ff A v (x , y ) dx dy 



Since the functions 9H/3pj satisfy the Lipschitz condition, we obtain 
from the first formula (1. 22), because of the inequalities (1. 27) and 
(1.28): 


’■ |£T ( 4 ; 


Av (x, y) d 


x dyj 


X , Y 


Y X 


a r 

+ Ti / / A v ( r, y) dy j dx -f T 3 / / Av (x, y) dx 


dy 


0 0 


0 0 


Cons equently , 


jr 1]i j < (Ti XY + T Z Y + T 3 X) u jAv (x, y) | 2 dx dy , 

G 

and where are definite positive constants. 

The functions 3 Z H./ 3p^3pp are bounded in the G region. Consequently, 


| r| | < (T 4 XY .+ T 5 Y + T 6 X) //i A v (x, y) ] 2 dx dy . 

G 

In this manner, the residual term in the formula (I. 21) satisfies the 
inequality 

|r)| — (Ames G + BX + CY) f I A v (x, y)| 2 dxdy , (1.29) 

G ' 1 

where A, B, and C are definite positive constants. If the function Av 
differs from zero on the circle of radius €, then from (1.29) it 
follows that 

| r| | ^ L,e J* J* Av 2 (x, y) dx dy , (1.30) 

g € 

where L does not depend on e . 
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4. Proof of Theorem 1. The Case of Linear Control System 


From the formula (1.21) for the increments of the functional 
and the estimate (1. 30) of the residual term within this formula, one 
can easily obtain the proof of Theorem 1. 

As a matter of fact, let for the sake of definiteness, v(x, y) be a 
control which is min-optimal according to S, and z(x, y) and u(x, y) - 
the solutions of boundary problems ( 1 . 1 ) - ( 1 . 2) and ( 1 . 5) - (1 . 6 } 
corresponding to this control. Then for an arbitrary permissible 
increment Av(x, y) the inequality AS - 0 is valid. Let us assume 
that there exists a point (£, r \ ) within the domain G in which the maximum 
condition is not satisfied, i. e. , there exists a control v 1 such that 

H (i, T|, p (£, r\), V 1 ) > H (£, T|, p (£, r|), v(4, T|)) . (1.31) 

Since the functions z(x, y) and u(x, y) are continuous and z x , Zy, 
and v(x, y) are sectionally continuous, there exists a closed region 
G 1 C G containing the point (£, r \ ) in which the left- and right-hand 
sides of the inequality (1. 31) are not continuous and, consequently, 
are likewise not uniformly continuous. If (£, r\) is the point of dis- 
continuity of the control v, then it may be obviously related to the 
boundary of the domain G 1 . It follows from the inequality (1. 31) that 
one may specify a number 8 > 0 for which 

H (x, y, p (x, y), v 1 ) - H (x, y, p (x, y), v(x, y)) >6 (1. 32) 

in all points (x, y) CG f C G 1 , where G € - circle of radius €. Let us 
introduce the control 


( v( x , y) at (x, y) C G e 
V 2 (X, y) = l 

( v 1 at (x, y) C G € 


Then because of the relationships (1.21), (1.30), and ( 1 . 32) 

as = - jf/V. y> y)’ vl ) - H ( x > y> p( x > y)) j dx dy - p 


//- 


Is 


< - I I 6dx dy + | P | — - J I j 6 - e L, | A v (x, y ) | 2 { dxdy , 
Gf G: 


where A v = v 1 - v(x, y). Since the function Av is bounded, one may 
choose the number € so small that the expression within the square 
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brackets in the first part of the last inequality may be positive. Then 
AS will become negative which contradicts the assumption about the 
min- optimality according to S of the control v(x, y). This proves the 
theorem. 

Formula (1. Zl) for the increments of the functional together with 
the formulas (1. ZZ) for the remainder term allow the establishment 
of a more general result for the linear boundary problem. 

As a matter of fact, let the control process be described by the 
boundary problem 

m 

z ixy = 2 | c ik( x ’ y) z kx + d ik( x ’ y) z ky + gik( x > y) z k| + fi(v) , 
k=i 

z i(o. y) = cpi(y)» z i( x > o) = 4 j i(x) , ^.( 0 ) = 4 ^( 0 ), (1.33) 

i = 1 m, 

and let search for the control over which the functional S attains its 
minimum (maximum) value. In such a case 


m 


H 


(x, y, p, v) = 23 u i I c ik z kx + d ikzky + gik z k| 

i, k=i 


m 


+ 23 u i fi > 

i= l 

and the function forms the solution of the boundary problem 


m 


u ixy = E 

gki u k " 

k= 1 

m 

U ix (x, Y) 



k-i 


m 

u iy (X, y) 



k=i 


d 

"dx 


_d_ 

dy 


( 1 . 34 ) 


Uj (X, Y) = - A i , i = 1 , . , . , m . 
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Since according to what was proved earlier 


Au! (x, Y) = Auj (X, y) = 0 
(see formulas ( 1 . 13)) we find 
Aui (x, y ) = 0 . 


F urthermor e , 

9H(x, y, p, v + Av) 8H(x, y, p, v) 


9wj 


0 , 


w - (zi , . . . , z m , . . . , ziy, . . . , z my) 

Consequently, rjj = 0. We now calculate q z . W e have 


3H (x, y, p, v) 
Owx Sw^ 

'This means that 


= 0 , i , k = I , . . . , 3m . 


e e // 

i = i k- 1 G 1 


0 2 H (x , y, p+OAp, v+Av) 

0Uj 0Wp 


3 2 h (x, y, p + GjAp, v + A v) 
8ui8w]^ 


A u-Awj^. dx dy 


Since Auj(x, y) h 0, it follows that r\ z ----- 0. Consequently, in the case 
under investigation, the formula (1.21) takes the form 



y> 


p, v + Av) - H (x, y, p, v) 


dx dy 


(l. 35) 


By means of the last formula one can easily prove: 

THEOREM 2. For the permissible control v(x, y) of the boundary 
problem (1.3 3) to be locally min-optimal (m ax -optimal ) according to S, 
it is necessary and sufficient that it satisfies the condition of maximum 
(minimum ) . 
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5. The Control of a System by Means of Boundary Conditions 

Until now we assumed that the control is carried out only by 
means of the function v entering into Equation (1. 1) or (1. 33). The 
boundary values (1. Z) of the function were considered fixed. How- 
ever, the method which is presented permits the solution of a more 
general problem. 

Let the control process be described by the system of equations 
(1.1) while the boundary values of the function z^ are specified not by 
the conditions (1. Z) but by means of the differential equations 


Ziy (0, y) = CD i (y. Z1 > • 

1 ^ 

. . , Z m , V ) , 


z ix ( x > °) - ( x > zi , . 

■ • . z m > v2 ) 

(1. 36) 

initial conditions 



Zi (0, 0) = z- , 

i = 1 , . . . , m , 

(1. 37) 


where the functions and 4^ are continuous over y and x and twice con- 
tinuously differentiable over the remaining arguments, vi and v^ are 
control parameters taking values from the domains V 1 and V 2 of s- and 
t-dimensional euclidian spaces, respectively. 

The presence of parameters within Equations (1. 36) permits the 
control of the process by means of boundary conditions. To the per- 
missible controls within Equations (1. 36) we relate also the sectionally 
continuous functions v 1 (y) and v 2 (x) with values in the regions V 1 and 
V respectively. It is known (see, for instance, work of Sansone , 
pp. 16 and 17) that each pair of permissible controls v 1 (y) and v 2 (x) 
determines by means of Equations (1.36) and conditions (1.37) a 
unique pair of absolutely continuous functions z(0, y) and z(x, 0). We 
will understand in all what follows under permissible control within the 
boundary problem (1. 1 )-(l . 3 6 ) - ( 1 . 37) the function 

w (x> y) = (v (x, y), V 1 (y), v 2 (x)) , 

whose components are sectionally continuous functions with values in 
the domains V, V , and V , respectively. Consequently, to each per- 
missible control co(x , y) there corresponds a unique solution of the 
boundary problem (1. 1)-(1. 3 6 ) - ( 1 . 37) with the same smoothne s s con- 
ditions which we introduced for the boundary problem (1. l)-(l.Z). 

Let us introduce the notation 


Z1 



“ ( z i ) • • • f » • • • i ^xxi ^ } 

m 

Hi (y, q, v») = E u m (y. z. v 1 ) , 

i=i 

m 

H 2 (x, q, v 2 ) = S u^i (x, z, v 2 ) 


The function ui is defined by means of Equation (1.5) and the additional 
conditions ( 1 . 6). 


In the general form one is not able yet to solve the optimum problem 
with boundary conditions (1. 3 6 ) - ( 1 . 37). However, it may be solved 
using the above described method if the following conditions 


9fk ( x > Y> z > z x > z y > v ) 

^17 

9fk ( x > y. z » z x * z y > v) 

^ z ix 


= 0 


y 


x= 0 


(x, z, v ) 
dz{ 

9cpk (y, z, v 1 ) 
9z^ 


k, i = 1 m , v c V , v 1 C V 1 , v 2 C V 2 


(1. 38) 


are satisfied. 


Thus in. what follows we assume that the conditions (1. 38) are 
satisfied and, consequently, that for an arbitrary function u(u x , . . . , i^) 
the equalities 


9H (x 

. y> p, v) 

_ 9H 2 

(x, q, v 2 ) 


9Ziy 

I 

o 

II 


9H (x 

. y. p. v) | 

9Hi 

(y. q. v 1 ) 

9z ix 

x= 0 

3z^ 


1 2 

are satisfied no matter what the values of v, v , and v from the 
1 2 

regions V, V , and V , re spectively are. 

We will say that the permissible control oo(x, y) in the boundary 
problem ( 1 . 1 )-( 1 . 36)-( 1 . 37) satisfies the maximum conditions if 
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(1. 39) 


H (x, y, p (x, y), v (x, y)) ((=)) sup H (x, y, p (x, y), v) , 

vCV 

Hi (y, q (0, y), v 1 (y)) ( = ) sup H] (y, q (0, y), v 1 ) , 

v*C V 1 

H 2 (x, q (x, 0), v 2 (x)) ( = ) sup H (x, q (x, 0), v) , 

v 2 CV 2 

where z(x, y) and u(x, y) are the solutions of the boundary problems 
(1. 1 ) -( 1 . 36 )-( 1 . 37) and ( 1 . 5 ) - ( 1 . 6 ) corre s ponding to the control to (x, y) - 
(v(x, y), v^y), v 2 (x)) while the symbol (-) indicates an equality which is 
valid almost everywhere within the domain of the change of the argu- 
ment. The conditions of minimum are defined in an analogous manner. 

THEOREM 3. For the permissible control co(x, y) in the boundary 
condition (1. 1)-(1. 36)-(l. 37) to be min-optimal (max-optimal) according 
to S, it is necessary that it satisfies the conditions of maximum 
(minimum) . 

The proof of this theorem is carried out following the same scheme 
as in the case of the proof of Theorem 1: one first finds a formula for 

the increment of the functional followed by the estimate of the residual 
term, and only then one proceeds to prove the theorem. 



For the establishment of the formula for the increment of the 
functional, we take an arbitrary permissible control co(x, y) and 
denote by z(x, y) and u(x, y) the solutions of the boundary problem 
(1. !)-(!. 36 ) - ( 1 . 37) and ( 1 . 5 ) - ( 1 . 6 ) corresponding to this control. Then 
the equality 



go = 


fflt 

G L i=i 


ixy " H ( x > y> P’ v ) 


dx dy 


X r m 
/ £ u i < x > 

J o L i = l 


0) z ix (x, 0) - H 2 (x, q (x, 0), v ) 


dx 


Y r m 

+ f I F Ui (0, y) zi y (0, y) - Hi (y, q (0, y), v 1 ) 


/* £ u i (°» 
J o L i = i 


dy 


0 


is valid. 
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Let us denote by Aco the arbitrary permissible increment of the 
control w(x, y) and by Az and Au the increments of the functional 
z(x, y) and u(x, y) corresponding to this control. It is clear that 

A I = I |p + Ap, go + Ago | - l|p, co| = 0 . 

To transform A I we start from the equation 


ffp^xy dy dx + f p (x, 0) q x (x, 0) dx + f p (0, y) q (0, y) dy 

G 4 *y» 

= ff qp * y dy dx - f q ( x > Y) Px < x > Y) dx 


1 

- f q (X, y) p y (X, y) dy + p (X, Y) q (X, Y) 

J o 

+ p (0, 0) q (0, 0) , 


(1.40) 


which is valid for arbitrary twice sectionally continuously differentiable 
functions p and q which equation may be derived from the Green's 
formula (1 . 16). 


After putting into Equation (1.40) p = A Uj and q = A zj and taking 
into account the conditions (1.38'), we find 



X m 

Zixydx dy + f 

o i = i 

Aui (0, y) A zi y (0, 


Aui (x, 


y) dy = 


0) A z^ x (x , 0) dx 




A zj 


. 3H 3H 

+ A ap: A2ix+ A aV" Az iy 


IX 


L y 


X 

m (r 

f 

E 


Au ix (x, Y) + A 


dx dy 

9H (x, Y, p (x, Y), v) j 




A zj (x, Y) 



9H (x, 0, p (x, 0), v) 

A ^ Azi (x, 0 )>dx 


9z 


iy 


Y, m 

- Z 

-i—l 


Au^y (X, y) + A 


9H (X, y, p (X, y), v) 


- A 


o i = H 

9H (0, y, p (0, y), v) 


9z 


IX 


Az i (X, y) 


9 zi x 


A z i (0, y) > dy 



m ' A 9H A A 9H A A 9H A 

A 9 z" AZi + A 0 Z . " A z ix + Az iy 

G i- i L 0Zl z ix oz iy y 


dx dy 


♦ X r t. A 8H » <«■ *,.<»■ tl A zi (x, 0, dx 


/ Z 

*0 1 = 1 

Y m 


9z; 


+ / 

J o i = 1 


9Hi (y, q (0, y), v 1 ) 


0 Z; 


A Zi (0, y) dy 


On the other hand 


X 


J/ S Au i Az ixy dx + 


/ m 

E i, L «, 0) A,. u 


(x , 0) dx 


i = l 


i=l 


y m m 

j Aiij (0 , y) Az. y (0, y) dy = jj ^ A Ai^dx dy 


X 


+ I 7 A 8H ^ (X ’ q (X ’ 0) ’ v2) Au i( x, 0) dx 


J 111 

Z 

0 1 = 1 


9 u[ 


f V A 8H - (y ' q(0 ' Y) ’ V ‘> A„ i( 0, y)dx . 

0 i^l 3U; 
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From the two last equations we obtain 


//£ 

M i= i 


3 m 

Az ix dx dy + I ^ Aqfx, 0) Az ix (x, 0) dx 
o i= l 


f 


4m 


+ I 53 Auj(0, y) Az iy (0, y) dy = — 


o 1=1 


1 r( 4 ™ 9H 
7 I J 13 A — — A p 

2 [ j Gi=i 9 Pi 1 


. dx dy 


X 


J a 9H 2 ( X > qf*’ °)> V 2 ) 


A. 


i= l 


9qi 


Aq- (x , 0) dx 


f] ‘f A 8H , (y ' Ih Vl) Aq i( 0, y) dy 

o mi 


(1.41) 


Furthermore, using the same method as during the derivation 
of formulas (). 18) and (1, 19), we find 


G i 


u • A z • dx dv + 
i lxy 7 


X 

/ m 

Z “i(‘. 


0) A zi x (x, 0) dx 


o i=l 


\ m m 

I X u i ( 0> y) Az i y (°> y) d y = - 53 A i Az i( x - Y ) 

< i = i i - i 


ffs [8H 9H 


A X 9H A 

A z- + Az- 

1X 9 z- iy 

ix az iy 


dx dy 


X 


I 


m 


o i = l 


9H 2 (x , q(x, 0), v d ) 


9zi 


A?q(x, 0) dx 


Y m 

/ E 

^ 11 


9H . (v, a 10, vl. v 1 ) 


9Z; 


A z^ (0, y) dy 


[ i. 
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X 


re m " m 

JJe -i 


z ixv dx d Y + I 2 A M*> 
'G i= 1 o i = i 


0) z ix (x, 0) dx 


/ m 

jC A Ui(0, y) Z iy (0, y) dy = 

n 1 — 1 


Xfe 5- *•' 

G 1=1 1 


dx dy 


X 


+ , f 0). ll> Aui(x , o) dx 

da, 


/ ill 

X 

1= 1 


I 

/ m 

£ 

n 1-1 


9H x ( y , q(0, y), v 1 ) 
3 u i 


A u i (0, y) dy 


(1.43) 


Taking into account Equations (1.41), (1.42), and (1.43) and the 
fact that AI = 0, we find using the same method as during the proof 
of Theorem 1 the formula for the increment of the functional 


AS 



y, p, v + Av) 


H ( x, y, 


P. 


v) 


dx dy 



q(x, 0), 


q(o> y), 


v 2 + A v 2 ) 


v 1 + A v 1 ) 


H 2 (x, q (x, 0), v 2 ) | dx 

(y, q ( 0 , y), v 1 ) j dy - q , 

(1.44) 


where 


h - h + h + h 

1 2 3 


h 

1 


— rf ( f y > P> v + A v) - 9H(x, y, p, v) 

2 i = 1 I i 9pi a Pi 
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+ £ 

k= l 


r 0 2 H (x, y, p + 0Ap, v 4 A v) 


3 Pi 3p k 


3 H(x, y, p + 9j Ap, v 4 A v) 


9 p i 9p k 


A Pi A Pk | dx d y > 


x 

J 


2 m 

£ 

i- l 


3H 2 (x, q (x, 0), v 2 4 A v 2 ) 


3H 2 (x, q(x, 0), v ) 

3qi 


A q i (x, 0) 


+ £ 

k “ l 


r 3 2 H (x, q ( x , 0) 4 9 Aq, v 2 4 A v 2 ) 


9q i 9q k 


3 2 I4 2 (x, q(x, 0) 4 9 2 Aq, v 2 4 A v 2 ) 


Y 


n = rr 
3 2 



3qi 3q k 

3H i ( y > q(°> y)> yl + A v x ) 

3q t 


Ac U Aq kj dx > 


3R 1 ( y , q, (0, y), v 1 ) 
3qi 


A qi(0, y) 


4 


2 m 

£ 


k= i 


3 2 H ! ( y , q ( 0 , y) 4 9 4 Aq, v 1 4 Av 1 ) 

9q L 9q k 


3 2 H 1 (y, q ( 0 , y) 4 9 5 Aq, y 1 4 Ay 1 ) 
3qi 3q k 


A qi A q k dy . (1.44') 


Let us now estimate the residual term r\ in the formula ( 1.44). 
The magnitudes r\ and rj/? / are defined by the values of the function 
z and u on the boundary of the domain G. From Equations ( 1. 36) and 
the conditions ( 1. 37) we obtain because of the Lipschitz conditions 


28 


y m 


22 | a ( o , y)| < n |*22 l Az i(°> y)l d y + p f 2 ^ l Av k*(y)l dy > 
i=1 i i=1 •• i k=i 

m * m r t 

22 I A Zi (x, 0)1 < NjJ 22l Az i( x ’ °)l dx + Pj J 22 l Av k 2 ( x ^l dx • 


From this, according to the lemma mentioned before, it follows that 

Y 


|Azi(0, y)| ^ M 0 J | Av. 1 (y)| dy , 


0 1=1 
X 


|Az i (x, 0)| < Mj /Si Av 2 (x)| dx 


0 1 ~ 1 


We introduce the notations 


m 

a(x, y) = 22 


A Z- 


IX | 


1 = 1 


A 


v 


k = i 


x)| dx 

• 

( 1-45) 

m 


m 

= Si 

Az iyl ■ 

y = 52 i Az ii > 

i-1 

i = l 

Z! a 

v k‘ (y) l’ 

l Av2 ! =52 l Av k 2(x 


Since the functions h satisfy the Lipschitz conditions, then like in the 
case of the derivation of inequalities ( 1. 23) , we obtain here 


y r| 

a (x, y) < N 2 f a (x , y) dy + N 2 J [q (x, y) + Y(x, y)| dy 
0 0 

Y X 

+ N 3 / Av(x, y) dy + N 4 J Av 2 (x) dx + N 5 Av 2 (x) , 
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x b 

p(x, y) < M 2 | 6 (x, y) dx + M 2 j fa(x, y) + Y (x, y) ] dx 

o *o 

x f Y r 

+ M 3 J A v (x, y) dx + M 4 / Av'(y) dy + M 5 Av'(y) , 

*o *o 

x y 

Y(x, y) < fa (x, y) dx , y (x, y) < f (3(x, y) dy , 


0 < x < £ < X , 0<y<r|<Y 


From this we find: 


a (x, y) < N 6 J |p(x, y) + Y (x, y) | dy + N ? f A v ( x, y) dy 
0 o 


+ N g I Av 2 (x) dx + N 9 Av 2 (x) , 


P(X, y) 5 M 


J |a (x, y) + y(x, y) ] dx + M 7 J A v(x, y)dx 


X 

J* Av'(y) dy + MgAv'fy) 


Taking into account the estimates for the function Y, we find . 


r 

a(x, y) < N iQ J (3 (x, y) dy + N ? J Av(x, y) dy 


+ N s J Av 2 (x) dx + N 9 Av 2 (x) , 
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r x r 

(3 (x, y) < M 10 In (x, y) dx + M 7 I A v (x, y) dx 


M, 


I Av*(y) dy + M 9 Av 1 (y) 


(*) 


From these inequalities we obtain 

i n 


X. Y 


a(f, n) - N n I* | a(x, y) dx dy + N [2 f f Av(x, y) dx dy 

*0 “o *b *0 


Y Y X 

+ N, J 4 «(|, y) dy + N 13 j Av‘(y) dy + N a J Av 2 (x) dx 


+ N,Av 2 (f) , 


Vc 

P(£, n) < M n II 6(x, y) dx dy + M 12 l J Av (x, y) dx dy 


0 o 


0 0 


x Y. 

+ M 7 I Av(x, r|) dx + M. 


r x r 

I Av 1 (y) dy + M J3 / Av 4 (x) dx 


+ Mg A V 1 ( T]) 


Integrating the first of these inequalities over £ between the limits 
of 0 and £ and applying the above mentioned lemma, we find 

L rr \ 

Y < fa(£, t|) d| < N x4 JJ Av(x, y) dx dy + N I Av‘(y) dy 
o G o 

X 

+ N f Av 2 (x) dx 
16 J 0 

In an analogous manner, from the second inequality we obtain 
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J p(s. n) dri < m 14 IJ Av(x, y) dx dy + M lg J Av (x) dx 
*o G o 


+ M u j' Av*(y) dy 


From this and the inequalities (*), we find 

cc Y 

| A z- (x , y ) | < N h J / Av(x, y) dx dy + N J5 f Av 1 


(y) dy + N I Av 2 (x) dx, 


I A z ix ( x ’ y)| - nJAvC, y) dx dy + N lg f Av(x, y) dy 


+ N I Av dx 4- N_ a I 

19 J 20 J 

0 0 


y 

/ Av X dy + N 9 Av 2 (x) 


I A z iy (x > y)| y)dxdy + M 18 / Av (x, y) dx 


1 A 

Av 1 dy + J' Av 2 dx + V^Av'jy) 


(1.46) 


In an analogous manner we get 


|Au^(x, y ) | < Wl zi J*J*Av (x, y) dx dy + M 22 ^ Av 1 


(y)dy + M / Av 2 ( x) dx . 


1. 47) 


If Avj(y) - Av 2 (x) = 0 , while Av. (x, y) £ 0 , then from (!. 46) 

and ( 1 • 47) we obtain the inequalities ( 1.27) and (1.29). A fter establishing 
this fact, we go over directly to the proof of the theorem. 

Let for the sake of definiteness the permissible control c o(x, y) = 
(v(x, y), v (y), v 2 ( x ) ) be min-optimal according to S. Then during an 



arbitrary A co the inequality AS ^ 0 is valid. Let us assume that the 
theorem is not true. Then within the closed domain G one can indicate 
either a surface domain Gj within which the first equality ( 1. 39) is not 
fulfilled, or a segment of a straight line located on the boundary of G 
over which one of the two last equalities (1. 39) is not satisfied. 

In the first of these cases one can find the permissible control 
v C V such that 

H(x, y, p(x, y), v) - H(x, y, p(x, y) , v) > 0 for (x, y ) C . 

Then there exists a 8 > 0 such that 

H (x, y, p(x, y), v) - H(x, y, p(x, y), v) > 8 

for (x, y) C G^ C G x , where G e - a circle of radius e located 
adjoining the boundary in the interior of the domain G 1# Putting 
Av 1 = A v 2 = 0 and repeating the same reasoning as the one carried 
during the proof of Theorem 1, we obtain by means of the estimates 
( 1 . 46) and ( 1 . 47) that A S < 0. However, this contradicts the condition 
and indicates that the first equality is satisfied under the conditions 
of maximum. 

Let us investigate the second case. For definiteness we assume 
that the last equation (]. 32) is not satisfied. There exists then a 
control 7 2 C V 2 and a segment 1 of the y = 0 boundary of the domain 
G such that 

H 2 (x, q(x, 0), v2) - II 2 (x, q(x, 0), v 2 ) > 0 
for x CL Consequently, one can specify a number 6 > 0 such that 

Hz (x, q(x, 0), V 2 ) - H 2 (x, q, (x, 0) v 2 ) > 6 
for x C 1 CL where 1^ is a segment of length e. Let us put 

Avi z Av. 1 = 0 
x 

and study the auxiliary control 

Zj 1 (x, y) = (v, v 1 , v 2 ) 
where 
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_2 

V 


v for x Cl. 


v 2 for x Cl, 


Then the residual term 4 in Equation ( 1. 44) coincides with 4^ (see 


Equation (1.44*)) where 


A 2 -2 2 
Av = v - v 


and consequently, Av 2 differs from zero only for x C 1 € . 


Since the function 3H 2 / 3q^ satisfies the Lifschitz conditions and 


3 H 2 /3qj_3qj c are bounded, then because of the estimates (1.45) and 
( 1 . 47) we obtain 


- Me J | 2 l Av k 2 ( x )l) dx > 


1 , 


where M is a constant independent on €. By means of this estimate 
we can easily establish that AS > 0 and this contradicts the assumption 
about the min- optimality according to S of the control co(x, y). 


Theorem 3 is thus fully proved. 


Let now the control process be described by the system of linear 
equations 


m 


z. 


ix y 


= 12 | C ik (X) y) z kx + d ik (x ’ y) z ky + §ik (x ’ y) z k| + f i( v > 


k= l 


i = 1 , . . . , m , 

with the additional conditions 


(1.48) 


m 


z iy(° 


. y) = 12 c ik(°’ y) z k + ^(v 1 ) - 


k= i 
m 


ix (x, 0) = Y, d ik (x > °) z k< x - °) + 

k = l 


*i(0, 0) - 


1 , . . . , m . 


(1. 50) 



The special choice of the coefficients in systems ( 1 . 49) follows 
from the requirements (1. 38). Like during the proof of Theorem 2, 
we find that in the case under investigation the residual term r\ in 
formula (1.45) is equal to zero and, consequently, 


AS = 


, y, p, v + Av) - H ( x, y, p, v ) j dx dy 


X 


J [ H 2 ( x > C l> v2 + Av 2 ) - H 2 (x, q, v 2 )| dx 


o 

Y 


' j [ H 1 (y> V 1 + Av 1 ) - Hj (y, q, v 1 ) dy . 


From this formula follows the validity of the following theorem. 

THEOREM 4. For a permissible control co(x, y) in the boundary 
problem ( 1. 48) - ( ] . 49) - ( L 50) to be locally min-optimal (max-optimal) 
according to the functional S = € A^Zj (X, Y) it is necessary and 

sufficient that it satisfies the conditions of maximum (minimum). 
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Section II. OTHER PROBLEMS OF OPTIMUM CONTROL 
FOR HYPERBOLIC SYSTEMS 


We investigate the same problem concerning the' minimization of 
the functional S - 2 A^z^ (X, Y) in which the control process is described 
by the boundary problem (1 . 1 )-( 1 . 3 6 ) - ( 1 . 37) where z? , i = 1 , . . . , m 
are given numbers. The permissible controls defined in Section I are 
omitted by the requirement that the respective numbers zj (X, Y) belong 
to a convex set D of the space of the variables , . . . , z m . In this 
manner, in the problem under investigation, the permissible controls 
transfer by means of Equations (1.1) and ( 1 . 36 ) the point (z§ , . . . , Zj° n ) 
into points of region D. In what follows, we will assume that the convex 
region D contains internal points and that it is closed. 

For the solution of the problem like in the paper by Rozopoer 16 we 
introduce the function 

A (/.) - (A, y.) - J^Ajzj 

and denote by D* the set of points z* £ D* at which 

A(z*) - min A (z) 
z £ D 

If the set D* is not empty, then 

A (z ' ) ^A (z), z*CD*, z£D, 
and, consequently, the functional 

s = E A i z i < x - Y > > 


defined over the solutions of the boundary problem (1. 1 ) — ( 1 . 36 ) - ( 1 . 37) 
cannot take values smaller than A(z*). If there exists a permissible 
control which transfers the point z° into an arbitrary point of the set 
D*, then such a control is a min-optimal according to S. In such a 
case, the problem reduces to the calculation of controls transferring 
z° into a given domain. Such a problem will not be investigated in 
what follows, i.e. , we assume that there are no permissible controls 
transferring z° into D*. 
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1. The N ecessary Optimality Conditions 


We will say that the permissible control gj(x, y) satisfies the 
maximum condition relative to the giyen function u(x, y) if the condi- 
tions (1. 39) are satisfied, where z(x, y) is the solution of the boundary 
problem (1. 1)-(1. 36) -( 1. 37). 

THEOREM 5. If co(x, y) is min-optimal control according to S and 
z(x, y) - the corresponding solution of the problem (1. 1)-(L 36)-(1.37), 
then there exists a vector function u(x, y) relative to which the control 
cj(x, y) satisfies the maximum condition. 

Let 

w(x, y) = (v(x, y), v^y), v 2 (x)) 

be a min-optimal control according to S and z(x, y) - its corresponding 
solution of the boundary problem (1. 1 )-( 1 . 3 6 ) - ( 1 . 37). We denote by 
D“ (D + ) that part of the domain D for which 

A (z) « SA iZi (X, Y), zeD _ (A(z) £ 2 A izi (X, Y), z£D+ ) . 

The general part of these closed convex regions is the plane 
m 

• £ Ai ( Zi - Zi (X, Y)) = 0 , 

i=i 

which contains the point z(X, Y). Since the control co(x , y) is the min- 
optimal according to S, there do not exist permissible controls trans- 
ferring the point z° into the domain D~. Noting this fact, we introduce 
the variation of the control assuming that all the permissible controls 
are sectionally continuous. 

We choose arbitrary points 

(Xj, y j ) , i, j 5= 0 (x 0 = 0, y 0 = 0) 

in the domain G and denote by G- the rectangle formed by the adjacent 
points (x v , y^) in which the corner of the lower left angle is represented 
by the point (x^, y j ) . We establish the square 

Iij, x i+1 - T « X < Xi+1 , yj + 1 - T < y < yj +1 , 
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where the number t is chosen so small that for the given set of points 
(xi, yj) these squares do not have common points*. 

We take an arbitrary sectionally continuous vector -functions 
a — (x, y), (3^(x) and Yj(y) defined for x, y£ [0, l] and taking values in 
the domains V, V 2 , and V 1 , respectively, of the changes of the control 
parameters v, v 2 , and v 1 . We introduce functions 


v b ( x , y, a ik ) 


v d (y. Yk) = 


'b < x > Pi> = 


v(x 

. y) 

for 

(x, 

Y)£.h-i 

, k-i 



yk - y\ 

(x, 

y) c ii_i 


a ik 

( T • 

T ) f ° r 

’ k-i 

1 

V 

(y) 

for 

yC 

bk - T > 

yk] 

Yk 

( y V y ) 

for 

yC 

[yk - T > 

yk) 

2 

V 

(x) 

for 

x C 

[xj - T, 

Xi) , 

Pi 

cn 

for 

x e 

[x, " T, 

Xi) . 


The functions 

< x > y> a ik> Pi> Yk) = ( v b ( x > y> Q ik)> v b X <y> Yk)> 

v b 2 ( x > 8 i)) 

will be called the varied control, and £2 will denote the totality of all 
possible varied controls corresponding to all possible squares I^j and 
all possible functions aik, (3^ , and y^ of the above mentioned type. We 
denote the solution of the boundary problem (1. 1 ) - ( 1 . 36)-(l. 37) corre- 
sponding to the control cj^C a by z(x, y, co^). Then the function 

Az(x, y, u) = z(x, y, u b ) - z(x, y, w) 


is the solution of the boundary problem 


Az ixy( x » y. w) = A 


0H 

3ui 


, (x, y) C G, 


( 2 . 1 ) 


*In the case when the number of points (x^, yj) is finite, there are 
no doubts concerning the existence of such a t. 
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0 PJ 

A z iy (0, y,w) = A-^— i- , yC [0, Y| , Azi (x, 0,u>) 

a 9H 2 /- In vl 

= A — , x £ 0, X , 

du^ 1 J 

Az£(0,0,cj)= 0; i = 1 , m. 


(2.2) 


(2.3) 


Since Equations (2.2) are ordinary differential equations, then, accord- 
ing to results of the paper by Rozopoer 16 it follows from (2. 2) and (2. 3) 
that: 


y m 

6zi (0, y, w) -J' ^ 
o k= i 


cpi (y. z(Q. y>u),u) 

3z k 


5z k (°» y» w ) d y 


+ 


k 


L R i|vj’Vj 

j= i 


^ < y < yk+i> 


where 


x 


6z i (x, 0, to) -J 23 ~ '^ X ’ Q,0J ^ ,U> ^ 6z k (x, 0, co)dx 


o k= i 

1. 

I 

j=i 


f, 

+ S Q i | x j ’ Pj] ’ x l < x < x l + 1’ 


6z • 



R i[ y j. Vj] = J' [<Ti(yi> 2(0. Yy vl )> Vj(y)) - <Pi(yj. z (°> yj, V 1 ), v 1 (y))]dy, 
0 

1 

Q i( x j> Pj) f [+i( x j> z ( x j > °> v2 )> Pj( x )) 

0 

- ^(Xj, z(xj, 0, V 2 ), v 2 (x))jdx. 
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It was shown in the same paper that 


1 m 


6 Zi (x,0,w) = ^ 2 A i s ( x - x j) Q s( x j’ w b>’ 


J = 1 S= 1 
k m 


(2.4) 


Sz^O, y, co) - EE B is(y>yj) R s (yj- 0 J b); 


J=i s = i 

where the matrices A^ g and B^ s do not depend on the choice of the func- 
tions 6. and Y • . 

1 i i 

From the relationships (2. 1) it follows that 

Az.(x, y , co) = Az*( x, 0, co) + Azj(0, y,co) 
x y 

+//E E ‘uv ki (2. S) 


0 0 S = 1 


+ E if (x, y)GG |)k - I z ( 


k' 


Here we introduce the following notation: 

x j yk 

J ijk [^bj = f f F i(x, y, w, aj k , v)c!y dx. 

x i“ T yk" T 


-IE fj 


s , q= i oo 

/ k x j y v 


3 w + eAw,v b ) 

3 w s 3 w q s 9 y 


EE// E— 4^'° ii, ' v) Aw sdydx , 


J - 1 V - l X j_ T y v ~ T s - 1 


where 


F i(x, y, w,a JV , v) = f{ 


(x, v. w .a * 

x j ■ x 

yv - y\\ 

\ y jv 1 


T / 


T~jj~ V x> y» w > v )- 
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In an analogous manner we find that 


Az ix^ X ’ y,t ^ = 


+ 


Az. (x, y,oo) = 

l y 


+ 

with (x, y)CG i, k 
where 



Az ix^ x ’ 0>w ) + 


/• 3m 

/£ 

0 S - l 


di[(x, y, w, v) 

9w„ 


Aw dy 
s 3 


S E iv | w b|+ E i> 

V - 1 


( 2 . 6 ) 


a /a \ C V"' 9 f-(x, y, w,v) 
Az iy ( 0 ,y ,»)+j ^ Aw g dx 


0 s = 1 


Z E ij K + e l 


j=i 


( 2 . 7 ) 


- I 


h k 


forx ; < x< x ; +1 - T > 


;x,y,w,a ? , p ,v)ciy for x ; + j -t x < x / + 1 , 


for v, < y < 


'k+i 


x, y, w,Qj k , v) dx for y k+ i - t y< y k+ r 


' 3 m 

/" V"' 3 2 f:(x,y, w + 9 Aw, vu) . . , 

7 A 1 5w s 8w q — Aw s A Vy 


o S , q = 1 


for x, < x < x 7 . - t , 

l 7 + l 


y P 3m _ 

^ y £) 8 F i( x,y.w,,, |P ,V) Awsdy 

i y p - T s — i 


3 w c 


for x, , - t ^ x < x, , 

7 +i 7 + i , 
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r 3 m ^ _ 

L .i V [ l fj(*.y.» + eA "'''b> Aw Aw d 
1 z Zj J 8 w s 3wq s 3 


x 


s , q = l o 


Fi=, 


for y k < y< y k+1 -t, 


L: + 


l J 3m 

E / E 

j=l x i" T 55 = 1 


3F i (x, y ,w,a Lij ,v) Aw dx 


J w. 


for Yk+i - T ^ y< Yk+i- 

According to what was proved earlier (see unequalities (1.46)) one 
can specify a positive number N such that 

|Awjjx, y)| N t , 


and consequently, 


|E l | Njt 2 I E i J $ N 2 t 2 |F i | « N 3 t 2 


and uniformly over x and y 


xl. • , x • ^ 

lim E - lim 2 - = lim 0 . 

T T T 

T— ►-O T-*- 0 T— ►* 0 


Introducing the substitutions 


= Xj - X, r|T - y v -y 


and going over to the limit we obtain 

l i 

R ijp [ x j’Yp’ w b| = lim = f f [ f i(*j> Yp.w(Xj, y p ), a jp(e,ri)) 

0 0 0 


- f i( x j» Yp’ w ( x j »Yp) v ) d ^ dr i- 


One can show that the totality of Equations (2. 5), (2. 6), and (2. 7) 
is solvable and for all (x, y) not located on the lattice x - Xj, and y = y y 
there exist limits 
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lim 

T — ► 0 


Az i( x -y- u) = 6z. 


lim - 

T— ► 0 


__ Sl . x> „ m A» ly (x.y.») = 6 


T — ► 0 


where 


a6z i c a8z i 

6z • =— — L. , 5z: v = — 


-ix ax 


Ji y 9y ' 


Dividing these equations by t going over to the limit for t — *-0, 
we find: 


6z^(x,y,co)=: 6z|(x,0,co) + &Z|(0,y,w) 

// E 3 fi, L y ; w ’ v) 6wsdx dy 


0 0 


S - 1 


l k 

+ E S Ri jpi x r y p'“i- 


j= i p= 1 


> (2-8) 


y 


r i ^ 8f 1 *(x,y,w,v) c j 

6w p dy ’ 


/ ^ 

b = i 


r 9 f:(x,y,w,v) c , 

6z iy (x,y,o) = 6z iy (0,y,») + f V 6w s dx 


S = 1 


for 


x z <x<x Mi = x > yk < y< yk+ 1 - Y - 

From the way the functions 6z^(x, 0, w) and Szj(0, y, u>) were defined 
it follows that 
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Sz^x, 0, w) - 6zj(0, y, co) - 0 


for 


0 ^ x < x x , 0 ^ y ^ yi, 

and, consequently, from (2.8) it follows that 
Sz.(x, y, cj) = 6z ix (x, y, w) = 6z iy (x, y, co) = 0 

for 


0 ^ x< x 1( 0 y < yj. 

Further, from the relationship (2.4) and (2.8) we obtain 

* Y 3m 


m 


6zj(x, y, u) ^ A is (x, Xj) Q s (x!w b ) + ff ^ — v l 6w g dy cl 


OX 


s 1 


Xl 0 S 1 


m 


6z ix (x, y,w) - ^ A is , (x, xj) Q s (x^) f J ^ ~ 6w s c] y 


S - 1 


0 s i 


3 m 


67. iy (x, y, u) f 


9fj(x,y,w,v) s dx 

8 Wo 


x 1 S - t 


for 


Xi < x < x 2 , 0 < y < yi . 

Solving this system, for instance, by the method of successive approxi 
mations, we find that the function 8z^ may be represented in the form 

m 

8z^(x, y, to) = ^ ] A\ s (x, y, xi) Q s ( x i ,^>b) , xx < x < x 2 , 0 ^ y ^ y x 

s -- 1 

where A*- s (x, y, xi) is the fully defined function not depending on the 
choice of the functions a^j, p., and yj. Continuing these deliberations, 
we define 



m l 

5 Zi (x, y, w) =2^ ^- s (x, y, x j) O s (xj, w b ), 

s = i j =1 

0 ^ y ^ yi, xj < x ^ xj+i = X. 

We find in an analogous way that 

m k 

6 Zi (x, y,w) = B f s (x* y> y p ) R s (y p > w b )- 

s = i p=i 

o^x^xj, y k < y ^ y k+1 = Y - 

From these relationships it follows in particular that 


( 2 . 9 ) 


(2. 10) 


m 

6 Zi (x; + 0, y,u) - 6 zj (xj - 0, y, w) = [A ? is ( x ? + y. x j) Q s ( x r W ^ 

S = 1 

- A /s 1 ( x z - °> y> Q s ( x z-i' w b)| > 

m 

5 z i(x, y k + 0,w) - 6 Zi (x, y k - 0,u>) = [b^x, yk + °> yk) R s (yk> w b) 

S- 1 

-Bb- J (x, y k - 0, y k _ 1 )R s (y k _ 1 , Wb)|. 

Consequently, functions oz^(x, y, co) which are defined by the formulas 
(2. 9) and (2. 10) are generally speaking discontinuous along the lines 
x = xi and y = y^,. 


Continuing analogous deliberations, we get 


m l k 

•i". -’-£££ l c ij ps (x ’ y ’ x j’ y p ) s (x j’y p > w b) 

S= 1 j- 1 prl 

+ D h < x ’ y> x j) Q s ( x j> w b) + F is P ( x > y» y P ) R s(y p . w b)| - 


x z < x < x z+i = x > y k < y ^ yk+i = Y • 

Substituting in this equality x = X, and y = Y, we obtain finally 
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( 2 . 11 ) 


m I k 

Sz^X, Y , w) = EEE | C ijps< x j« V S (x j y p ,W b ) 

s=i j=i p=i 

+ D^ ig (Xj) Q g (xj,u b )+ F? s (y p ) R s (y p ,wb)| , 

where the constants C- . , D. 1 ^ , and F? do not depend on the choice 

Ijps IS IS r 

Of ay, (3 i , and Yj- 
The point 


z( X, Y) + 6z(X, Y, (o) 


corresponding to the arbitrary variation co^ £ ft of the control to goes 
over a certain set II within the space of the variables zi, . . . , z m . In 
the same manner as it was done in the paper by Rozopoer 16 one can show 
that it is convex and that its arbitrary internal point cannot belong to 
the internal part of the set D . It follows from this that through the 
point z(X, Y) one can draw the plane 


m 

^ a i {z i - Zi (X, Y)) = 0, 

i - 1 


( 2 . 12 ) 


dividing the sets II and D where the signs of the coefficients a.^ may be 
chosen in such a way the II is in the half- space 

m 

^ a.(zj - Zj(X, Y)) 3= 0. 
i= i 


Consequently, for an arbitrary co^ 
y; a i 6z i (X, Y,w) 550, 

i. e. , 

£ a i Az i (X, Y , to) 

lim * ^ 0. 

t — +-0 T 

We introduce auxiliary functions m by means of Equations (1.5) 
and additional conditons 
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( 2 . 13 ) 


u ix^ x ’ Y ) = 


3H(x, Y, p(x, Y), v) 
9z iy 


U iy( X ’ 7) 


9H (X, y, p(X, y), v) 


3Z; 


IX 


u.(X, Y) = -aj. 


Using the same method which we applied above, we can obtain a formula 
for the increment of the functional 

S = Ea i z.(X, Y) 

in the form (1.44) and, consequently, the same method may be used to 
show that the conditions of the maximum ( 1 . 39) are necessary in order 
that the permissible control co(x, y) realizes the minimum of the functional 
S. However, S attains its minimum over the min- optimal control 
according to S. 

Theorem 5 is thus fully proved. 

It is obvious that the statement just proved remains valid even in 
the case when the control process is described by the boundary problem 
( 1 . 1 )-( 1 . 2 ). 


If the control process is described by a linear boundary problem 
(1.48)-(1.49)-(L 50), then one encounters as valid 

THEOREM 6. Let z(x, y) be the solution of the boundary problem 
{ 1. 48)-( 1 . 49)-( 1.50), corresponding to the control o>(x, y) and satisfying 
the condition z(X, Y) = z 1 . Then, if w(x, y) satisfies the condition of 
maximum (minimum) relative to the functions u^x, y) taking the boun- 
dary values 


Ui (X, Y) = -\A t - (J-B i (z 1 ), ^ 0, \>0, 

where B^z 1 ) a re the coordinates of the normal, perpendicular to the D 
hyperplane, then the control o(x, y) is min- optimal according to the 
functional 

m 

s A iZi (X, Y). 

i =i 

The proof of this theorem agrees almost completely with the proof 
of the corresponding theorem (see Theorem 4 in the paper by Rozopoer 16 
for ordinary differential equations. 
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2. The Use of Theorem 5 for the Solution of Certain Specific Problems 


The result just obtained does not, generally speaking, present 
a method for the establishment of the vector u(x,y). However, in num- 
erous particular cases this problem may be solved. Let us study some 
of these. 


1) The point z(X, Y) is located within the domain D. Then 
= Aj_, since an arbitrary plane in addition to the plane 
(Z. 1Z) crossing the point z(X, Y) cuts the region D and, 
consequently, cannot separate D and II. 

Z) The point z(X, Y) belongs to the boundary of the domain D 
which is specified by the inequality F(z)^ 0. Then the 
boundary is specified by the equation F(z) = 0. If the 
function F(z) is differentiable, then the equation of the 
tangential plane through the point z(X, Y) has the form: 


m 

B^zi - Zi(X, Y)) -- 0, 

i= i 


Bi 



z =z(X, Y) 


Since the plane 


Sajlz. - z.(X, Y)) = 0 


likewise crosses the point z(X, Y), then 

a i = XA- + pBp 

where, without a loss of generality, one can assume that 
X > 0, p $ 0 (\ l + p 2 * 0). 

Since a^ is determined with accuracy up to a constant multiplier, then 
only one of the quantities X and p is independent. Since, according to 
the conditions of (Z. 13), u^(X, Y) - a. and F(z(X, Y)) = 0, we obtain 
m + 1 relationships 


u t (X, Y) - -XAj - p B l> F(z( X, Y )) 


i = 1 , 


m, 


(2. 14) 


for the determination of u^(X, Y) and of the quantities \ or p. Adding 
to (Z. 14) the conditions (1. 37), we obtain Zm boundary conditions for 


the Zm functions z\> 


\yi th F! 


J m> 


ui , 


4 m' 


These conditions, together 


(1.1), (1.5), (1.36), ( 1.39), and (Z . 1 3) , form a 
complete 1 ' system of relationships for the' d etermination of the optimum 
control and the corresponding vector functions z(x, y) and u(x,y). 
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Suppose that we are required to determine the minimum of the 
functional 

X Y 


i = j" J M x > y. z > z x> z y > v ) d y dx 


0 0 


under conditions that the function z(x, y) is the solution of the boundary 
problem (1. 1)-(1. 2), while the point z(X, Y) belongs to a certain convex 
domain D of the space of variables z\, . . . , z m . By introducing an 
auxiliary function z 0 through the relationships (1.4), we reduce the 
problem to the calculation of the minimum. . . [Apparently one line is 
missing, Note of the Translator^ . . . enters into the cylinder having a 
generatrix parallel to the axis z 0 . Since the variable zo does not enter 
into the right-hand side of Equations (1.1) and (1,4), we find that Bo - 0 
in the relationships (2. 14). For the functional under inve s tigation, 


Ai - . . . 




= 0, Aq - 1, which means that it follows from (1.5) and 


(1.6) that u 0 (x,y) = -1. In this manner, differential equations and 
boundary conditions of the problem under investigation take the form 
of the r elationships 


z ■ 


lxy 


9H 

9u. 


z i(o> y) = a> j_( y ) » z.(x, o) = ^(x), 


u 


_ 3H 
* x y 3z i 


d / 3 H \ d / 3H \ 
dx v 9z ix ) d y \ 9z iy ) ’ Uix X ’ 


3H 


3z • 


iy 


y = Y 


u i ( x ,y) 


3H 

3z. 


u. ( X, Y) = 0, H 


ixl 


x = X 


E v, - 


fo, 


1 - 1 


from which the auxiliary Equations (1,4) are excluded. 


3. Generalization to the Case of an Arbitrary Number of Independent Variables 

The formula for the increment of the functional S and the 
ensuring consequences may be generalized for the case when the control 
process is described by the Goursat problem with an arbitrary number 
of independent variables 18 . Nevertheless, to avoid a cluttering of the 
formulas with irrelevant details, we assume that the number of inde- 
pendent variables is equal to three. 

Thus let the functions 
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z^x), x = (X1.X2.X3), i = 1 , .... m, 

be specified through the relations 

9 3 Zi / 9 Zl 9z m 

8x18x28x3 0:1 9 xj 8x3 

i = 1 , . . . , m, 0 ^ < X k> k = 1,2, 3, 

and the additional conditions 

2 j( 0 , x 2 , x 3 ) = cp. 1 (x 2 , X 3 ), Zj(xi, 0, X 3 ) = (J)^ (xi, 3 

z.(xi, x 2 , 0) = CDi 3 (Xl, x 2 ), 

where the functions T contain mixed derivatives of the zj variables of 
an order not exceeding two. These functions are twice continuously 
differentiable over the totality of all arguments. The contr ol^par a- 
meter obeys the same conditions as before. The functions rpi are 
twice, sectionally, continuously, differentiable over their arguments 
and satisfy the natural matching conditions. Like in the previous cases 
we will assume that each permissible control has an associated class 
of functions within which the Goursat boundary problem can be uniquely 
solved. 


9 2 z 


m 


8 x 2 8 x 3 ’ 


I 


(2- 15! 


< 2 . 16 


As the criterion of optimality we chose the functional 
m 

s = A iZi (Xx. X 2 .X 3 ), 
i “1 


We introduce auxiliary variables ui and the function 


m 

H(x, w, v) = u^, 

i-i 


(2. 17 


where 



u nV z i , 


8 z i ^ z m ^ z, m \ 

Zm * 9 xj ’ * 9 x 3 * ’ 8x28x3 / 


is a vector with a number of components equal to N. The function 
m(x) is defined by means of the equations 



3 3 ui 

9H 

-E- 

k =1 

9 / 9H\ 

D 

— V- 

9 2 

/ 9H \ 

3xx3x23x3 3 

8xk \ 3z ix k / 

2 2-U 

j * k 

3xj3x^. 

\ 9z / 
\ ix jxk / 

1 ? 

i = 1, - 

. * , m, 






\ - • 

the auxi 

liary conditions 






3 2 ui 

9H 

9 / 

8H \x. 

9H 

/ 9H 

| for 

X 3 - x 3 , 

9xi9x 2 

9z ix 2 

9xj l 9z i f 

\ 1 x 1 x 2 / 

9x 2 

\ 9z ix 2 x 3/ 

9 2 Uf 

9H 

.+ 9 1 

i 3H \ + 

9 / 

' 9H \ 

for 

x 2 =X 2 , (2 

3xi9x3 

9z ix 2 

9xj 1 

l 3z ix,xj 

9x 3 \ 

^ 9z ix 2 x 3 / 

3 2 iii 

9H 

+ 9 1 

1 3H \ 

9 

/ 9H ' 

| for 

XJ = Xj , 

3x 2 3x3 

3z- 

1 x 1 

9x 2 | 

\ dz ixixJ 

8 x 3 

\ 9z ixix 3 y 


3u^ 3H 

~ 9 z ix 3 x 2 

9u t 9H 
9x2 9zj XlX3 

9uj _ 9H 

9 x3 9z ix\X 2 


for x 2 - X 2 , 


for X 3 = X 3 , 


for xj = Xj , 


x 3 = X 3 ; 

Xi = Xj, 
x 2 - X 2 , 


u^(Xi , X 2 , X 3 ) — — , i 1 » • » • » m. 


( 2 . 20 ) 


( 2 . 21 ) 


Equations (2.20) are ordinary differential equations. Consequently, 
for each permissible control, they determine together with the conditions 
of (2.21) uniquely the functions u^( xj , X 2 , X 3 ) , u 3 (Xj , x 2 , X 3 ) , and u^(Xi , X 2 , 
x 3 ). We now solve Equations (2. 19) with the additional conditions 


U^Xi, x 2 , 

x 3 )| 

Xl 

= X! 

U^Xl , x 2 , 

x 3 )| 

*2 

- X 2 

Ui(xi, X 2; 

. x 3 )| 

Xl 

II 

u^xi, X 2; 

. x 3 )| 

X 3 

= X 3 

Ui(X 1( x 2: 

. x 3 )| 

x 2 

II 

X 

rx; 

ui(Xi, x 2 

. x 3 )| 

X3 

= x 3 


u i (X 1 ,x 2 , X 3 ), j 
u i (x 1 , X 2 , Xj), I 
u i (X 1 ,X 2 ,x 3 ), | 
u i (x 1 ) X 2 , X 3 ) , S 

Uj(Xi, X 2 , x 3 ), | 
Uj(Xj , x 2 , X 3 ) , ) 


for x 3 - X 3 ; 


for x 2 = X 2 ; 


for xj = Xi. 
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Because of the assumptions made above, the functions u^ ( Xi , x 2 ,X 3 ), 
u[ (x 1} X 2 , x 3 ), and ui(X 1 ,x^,x 3 ) are uniquely determined. In this manner, 
the problem is in the last resort reduced to the Goursat problem: to 

find the solution of the system of Equations (2. 18) within the domain 
0< Xk^ Xk, satisfying the boundary conditions 


Ui(x 1; x 2 , x 3 ) |x 2 = X 2 = u 3 ( x i> X 2> x 3 ) , 

14 ( x i> x 2> X 3 ) jx 3 = X 3 = U|{xi, x 2 , X 3 ) , 

Ui(x!, x 2) x 3 ) |x t =Xi =ui(X!, x 2> Xj) , i = 1, 2 m. 


( 2 . 22 ) 


At the same time, one must keep in mind that the functions 
3H/ 8 w^ within the Equations (2. 18) and (2. 19) are d iff e r entiated over 
Xl , x 2 , and x 3 . This means that if one assumes that the class of per- 
missible controls consists of s ectionally -continuous functions, the 
following conditions must necessarily be satisfied: the right-hand 

side of these equations should not depend on the derivative functions 
v and on z XlXj , z x 2 x z > and z x 3 x 3 * if the right-hand side of Equations 
( 2 . 18) and ( 2 . 19 ) do depend on these quantities, then one must choose 
the functions v(x) having sectionally continuous derivatives for the 
class of permissible controls. 


Assuming that these conditions are satisfied, one can obtain 
the formula for the increment of the functional 


XiX 2 X 3 


AS = 


-Iff 


H(x, w, v + Av) -H ( x, w, v) dx 3 dx 2 dxi — rj , 


using the same method as presented above, with 


r\ = r]i + q z , 




N X!X 2 X 3 

'{ZSSS I 


3 H(x, w, v +Av) 3 H(x 


i - 1 o oo 


, w, v) 1 

3wi J 


Awidx 3 dx 2 dxi , 


N X t X 2 X 3 

iE Iff [ 

i is- — i r\ r\ ' 


1,K = 1 0 0 0 


8 Z H (x, w -f 9 iAw, v+Av) 

8 w[8wk 


l8 2 H(x, w +0 2 Aw, v-fAv) 
8 w^8wp 


A w{ Awk dx 3 dx 2 dxi 



From this formula one can derive the optimality conditions 
which may be formulated in the form of Theorems 1 and 2. If one 
assumes that the control is carried out by means of boundary con- 
ditions, one can obtain results which are analogous to Theorems 
3 and 4. 


4. The Control of a Process by Means of "Concentrated Controls" 

We assumed, in all the problems investigated above, that 
all the components of the vector v(x, y) are functions of two variables: 
x and y. However, the proposed method permits the solution of a 
problem where all the permissible controls v(x, y) may be presented 
in the form 

v(x, y) = (v 1 (x), v 2 (x, y), v3(y)) 

(some components of this vector are functions of only a single inde- 
pendent variable x or y). 

For definiteness, we investigate the problem of the minimization 
of the functional ( 1 . 3) when the process is described by the boundary 
condition (1. 1 ) -( 1 . 2). The formula for the increment of the functional 
( 1. 21) remains valid even in this case. Valid is also the estimate 
(1. 29) of the residual term of this formula. Consequently, the same 
method proves also 

THEOREM lb For the permissible control v(x, y) -( v 1 (x), 
v 2 (x, y), v3 ( y )) in the boundary problem ( 1 . 1 ) -( 1 . 2) to be min -optimal 
according to the functional ( 1 . 3), it is necessary that the condition 



H (x, y, p(x, y), v(x, y) + Av) - H(x, y, p(x, y), 


v(x, y) dx dy^S 0 


be satisfied for an arbitrary permissible increment Av, where p(x, y) 
is a vector corresponding to the control v(x, y) and which is fixed by 
Equations (1.1), (1.5), and the additional conditions (1.2) and (1.6). 


If in particular the permissible control depends only on a single 
variable (for instance, on x) and in Equations (1.1) we have 

f i ( X) y, z, z X) z y , v) =fi(x, y, z, z x , z y ) +f/ (x, v), 

then the conditions (2. 23) take the form 
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m 

fJT. ut(x, y ) [ fi(x, v ( x) +Av) - 4 (x, v)] dx dy ^ 0 . 


x = 1 


By introducing the notation 

m i J 

H 1 (x, u (x), v) f i( x > v )/ u t ( x > y) dy , 

i = i o 

we obtain the optimality condition in the following form. 

THEOREM r*. For the permissible control v(x) of the boundary 
problem (1. 1)-(1.2) to be min-optimal (among controls depending only 
on x) according to the functional (1. 3), it is necessary that 

H 1 (x, u (x), v (x)) (=) sup H 1 (x, u ( x ) , v) , 

vCV 

where the sym bol ( -) indicates equality valid for almost all x 1 s of 
the segment 0^ x -- X. 
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Section III. THE VARIATIONAL CALCULUS AND THE PROBLEMS 
OF OPTIMUM CONTROL 


The problems investigated in the present paper represent in essence 
problems of variational calculus. However, the classical methods can 
not be applied here since the control parameters may, in general, take 
values from a closed domain. In the case when the region of variation 
of the control parameters is open, one obtains from the principle of 
maximum the necessary conditions of the classical variational calculus 
for functionals with partial derivatives. 

Let us search for the minimum of the functional 


I - 


X Y 

If 


f(x, y, z, 


J x ’ 


z v » v ) d y dx > 


0 0 

which is defined over the functions z = 
relationships 


(z 


i> 


z m ) specified by the 


z ixy ( x > y) = v i» v - ( v i> • • • > v m)> z i(°> y) - Ti (y)> 

Z{ (x, 0) = vJjj (x), i = 1, . . . , m, 

where the control parameter v is chosen from the class of all sectionally 
continuous vector functions. 


Under the optimum control we understand the permissible controls 
found within the immediate vicinity of the function z (x, y) corresponding 
to such a. control. It is obvious that such a definition of the optimum 
control is a special case of the optimum control in the previous sense. 
Consequently, the principle of maximum remains valid and every opti- 
mum solution is also an extremal solution. The opposite is also valid: 
each extremal constitutes an optimum solution. 

For the calculation of such a solution we introduce an auxiliary 
variable z 0 : 

z 0xy = f ( X > y> Z ’ Z X> z y > v), z 0 (x, °) = z o(°> y) = °> 

and establish the function H: 


H = u 0 f + 2 UpV p . 
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Then the auxiliary functions uj. (x, y) are defined by means of the 
boundary problem 


af 

11 • _______ 1 T 

d 

/ 0f 

lx y 0 Zi 0 " 

dx 

\ 9z ix 

u ix ( x> Y ) 


0f 


[ 9z iy U 

u i (X , Y) - 

0 , 

i = 1 


u, 


’) -37(^17 


U 0 


y = Y 


» u i y (X, y) 


3 f 

37— u ° 

° Z 1X 


M3. 1) 


Jx = X 


- - 1 . 


From this we find that 
H - ^ 1 u p v p “ 


Since the function H reaches its maximum on the optimum control 
v (x, y) , we have 


e). 


v — v (x , y) 
Consequently, 
d 2 


(“* 


u; 


1X Y dx dy 


(- 11 - \ 

\3 z ixy / 


1L \ 

0V ^ K = v(x 


y) 


0. 


and because of Equation (3. 1), we find that the solution z (x, y) of the 
optimum problem under consideration satisfies the system of equations 
by Ostrogradskiy-Euler (see, for instance, paper by Akhiyezer 32 
p. 1Z2): 


9f d , 

f 9f \ 

d 

/ 0f ^ 

j + d2 j 

f 0£ ) 

dz^ dx ' 

l 9z ix/ 

' dy 

\ 9z iyy 

! dx dy ' 

\3 z ixy/ 


As assumed earlier, the function f has a second continuous deriva- 
tive with respect to the variables v lf . . . , v m . Since the control 
v (x, y) realizes the maximum of the function H, the quadratic form 
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a z f , x 

9v i 9v k 1 k 


E 


9 2 H 

0 v i 9 v k 


x i x k 


m 

E 

i, k = l 


is nonpositive. Consequently, from the condition of the maximum 
(1.7), it follows that everywhere within the domain G(0 - x - X, 0 ^ 
y < Y), except, perhaps, in points located on a finite number of lines 
with zero area, the inequality (Legendre condition) 


£ 


i,k = i 


a*f(x, y, z ^zx. Sy- Zxyj ^ fe0 , £ X- + 0 


9vi9v k 


1 = 1 


(3.2) 


representing the necessary condition for the function z (x, y) to be 
extremal minimizing the functional I, is satisfied. 

In the case when the domain of variation of the control parameter 
is closed, the derivatives 3H/9v^ can not become zero along the optimum 
trajectory z (x, y) and, consequently, the condition (3. 2) may even not 
be satisfied. As a confirmation of what was just said we investigate 
the simplest case. 

Let the control process be described by the boundary problem 

z X y = v 2 , z (x, 0) = z (0, y) = 0, 0 — x , y - 1 , 

where v - the control parameter, jvj - 1. As a criterion of optimality 
we use the functional 


l i 

If 

o o 


j xy 


dx dy 


(1,1) (f (x, y , 


'x’ ^y* 


v) = - v 2 ). 


It is easy to show that the min-optimal control according to S is 
v (x, y) = 1 , and consequently, during this control 



0, 


and the condition (3. Z) is not satisfied. 
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Section IV. OPTIMUM PROCESSES IN SYSTEMS WHOSE BEHAVIOR 
IS DESCRIBED BY PARABOLIC EQUATIONS 


1. The Formulation of the Problem. The Maximum Principle 

Let E n be a euclidian space of vectors x = (x,, . . . , x n ) and 

n ^ 

G - a bounded region within E with a boundary F belonging to the class 
a(*> (see work of Miranda 33 , p. 10), and (x) are the direction cosines 
of the external normal on the boundary F. 

Let, furthermore, an elliptical operator L “ (L x , . . . , L m ) be 
defined within the domain G by the formula 


m n 

E E 


ip 


3 2 y P 


Jly Z-f Z— f jk 3x k 

p = 1 j. k = 1 

ip 


(4. 1) 


where the functions aj k (xj, . . . , x n ) within the domain G + r belong to 
the class C^ 2 ). We denote by M = (M lt . . . , M m ) the operator defined 
by the formula 


t 

P = i j. k = l j = l 


, m, 


where 


Pi _ 


n a a P x 
3a jk 


E 


k = 1 


7Xv 


One can check directly that the equation 


m 


£/ 


m n 


i G 


( Z ,L iy - Vi M iZ )dx = ^ ^ j ( al g-t 


i,p = l j = l r^k=i 


1 



is valid. In the same manner, used before for an elliptical type equation, 
we can transform this formula into the form 


m t rn - 

£ J (ziLiy - yiMiz)clx = 2^ J ( z i p iy - yiQi z ) do- , (4.2) 


i =1 G 


i = i r 


where 


m 

p iy =E 

p = i 


,'P^lE- + b 


dl 


l P 


ipypj . Qi 


m 

z =E 

p=i 


• dz D 

— + d^pZp 


d\ 


l P 


(4.3) 


In formulas (4. 3) the directions lip are chosen arbitrarily pro- 
vided cos (n, lip) > 0 (n - external normal to F) and their direction 
cosines belong to the class on T. The directions \ip are chosen 

depending on the lip. 

Let us assume that the coefficients within the operator L, depend, 
in addition, on the variable t, 0 ^ t ^ T, and we investigate the control 
systems whose behavior is described by a system of equations of the 
parabolic type 

L t y = f(t, x, y, y x , u), 0 « t T, xC G (Li t y = -If 1 - p iy) > (4.4) 


where the function f = ( fi , . . . , f m ) is continuous with respect to t 
and is twice continuously differentiable over the other arguments, while 
the parameter u takes the values within a certain convex (open or 
closed) region U of the p-dimentional euclidian space. 


We assume further that the function y(t, x) = (y lf . . . , y m ) is 
defined by the system of equations (4. 4) satisfying, in addition, the 
following conditions 


Pi (t, x)y 


= <T i (t, x, y, v), xC r, 0 < t < T, 
y (0, x) = a ( x), xC G, 


(4.5) 


where the operators P[ are defined by the formulas (4. 3) in which the 
functions aj s (t, x), b{p (t, x), and a(x) are continuous, ^ satisfy the 
same conditions as fp while the parameter v takes its value from a 
convex (open or closed) region V of a q-dimensional euclidian space. 

The function oo(t, x) = (u(t, x), v(t, x)) will be called the per- 
missible control if all its components are sectionally continuous and 
u(t, x) and v(t, x) take values from the domains U and V, respectively. 
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In addition, we will assume that the discontinuity surface of the per- 
missible control [Translator's note: one line seems to be missing.] 

. . . or in the vicinity of its arbitrary point one can introduce a non- 
degenerate coordinate transformation 

t ( 1> ^)> i ~ 

since the discontinuity surface goes over into a portion of the plane 

ln= 0- 


If the discontinuities of a certain permissible control satisfy 
the first condition, then the boundary probalem (4. 4) -(4. 5) corresponding 
to this control splits into several such problems located, however, in 
regions adjoining one another along the discontinuity surfaces of the 
control. In such a case the problem (4. 4) -(4. 5) has a unique continuous 
solution ( see, for instance, work by Zagoriskiy 34 ), and this solution is 
not subjected to any further additional smoothness condition over the 
discontinuity surfaces of the control. 

If these surfaces satisfy the second condition, then one views 
as the solution of the problem (4.4)-(4.5) the vector function y(t, x) 
satisfying the system of equations (4. 4), the conditions (4. 5), and 
certain smoothness conditions over the discontinuity surfaces of the 
control. In its most general form, this problem apparently was never 
studied although its particular cases were investigated in numerous 
papers 35 ’ 36 ’ 37 ’ 38 ’ 39, 40 which supply various existence and uniqueness 
theorems concerning the solutions. Consequently, we will assume 
everywhere in what follows that the given functions in (4. 4) and 
conditions (4. 5), in addition to the above listed properties, satisfy 
also the conditions under which to each permissible control there 
corresponds a unique solution of the problem (4. 4) -(4. 5). 

Let u>(t, x) be a certain permissible control and y(t, x) the 
corresponding solution of the problem (4. 4) -(4. 5), and let be given 
the functional 


s =T / Q i ( x )yi (T,x)dx + J J Pi (t, x)yi(t, x)dxdt 
1 = i LG o G 

T ■] 

Yi (t. x )Yi (t. x )do- dt 

o r 

where a^, p^, and Yi are given continuous functions. 


(4. 6) 


60 



L,et us formulate the problem: among all the permissible controls 

one should find a control oo(t, x) (if it exists) such that the corresponding 
solution of the problem (4. 4) -(4. 5) realizes the minimum of the 
F unctional S. 

The permissible control go (t, x) over which the functional S attains 
its maximum (minimum) value will be called the max-optimal (min- 
optimal) according to S. The functionals of a more general type will be 
studied at the end of the paragraph. 

As it was mentioned above, the problem of optimum control pro- 
cesses described by parabolic equations are of definite theoretical 
and practical interest. Numerous papers 5,6 ’ 11 investigated certain 
problems for which the control is materialized by means of initial 
and boundary conditions and for the optimality criterion one chooses 
either the speed or the functional of the type 

i T 

I = j* [u( T, x) - u 0 (x) | 2 dx + y J* p z (t) dt, 
o o 


where u 0 (x) is a given function from L, z (0, 1), p(t) - control, and y - 
a non-negative constant. We already investigated this problem for the 
case when the control of the process could be materialized simultaneously 
by means of controls entering into the equations as well as into the 
boundary conditions. It is clear that the functional 


Si = 



Yi(t, x ) Yi(t, x ) + 51 

k = i 


ayi 9yi 

a ik ( t, x) + Pi ( t, x) dx dt, 


9x k 


9t 


with a(k and Pi - continuously differentiable function can be reduced to 
the form (4. 6 ). 


To formulate the condition of optimality we introduce the auxiliary 
function z(t, x) = ( zp . . . , z m ) by means of the boundary problem M adjoint n 
to ( 4. 4) -( 4. 5): 


m 

M it z = -E 

S = 1 


9fs (t, x, y, y x , u) 
9y t 


n d /9f s (t, x, y, y x , u) 


k = 1 dxv 


ay,- 


ix k 


+ Pi(t, x), XC G, 
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(4.8) 


m 

Qi(t, x)z = 53 
s = 1 


acp s ( t ’ x > y> v ) n af s ( t > x > y> y x . v) 

+ 53 X k (x) | z s 

3y i k = i 8y ixk 


where 


] 


- Yi ( t, x), xCT, z i (T, x) = — Qi(x), xC G, i = l,...,m, 


A>r 3 z i x , 

M it z = ^ + MjZ, 

11 at 1 


are operators defined by the relationships (4. 3), the functions 04 , 
(3j, and yi enter into the definition of the functional S, and Xp(x) - 
the directional cosines of the external normal to G at the boundary F. 
For the boundary problem (4. 7) -(4. 8 ) to be solvable it is necessary 
that the functions (4 and y^ be connected with the matching conditions. 
In what follows we assume that these conditions are fulfilled. 


We introduce the notation 

/ 3yi 3y m \ 

w = (z L> . . . , z m , Yi. . . . . y m , — , . . . , , 

m 

P = ( z i > • • • > z m> yi. • • • . Ym)> H (t, X, w, u) =53 z{ fj (t, x, y, y x , u), 

i = i 
m 

h(t, x, p, v) =53 z i Cpi ( t, x, y, v). 
i = l 


Then the boundary problems (4. 4) -(4. 5) and (4. 7)-(4. 8 ) may be written 
in the form 


9 H(t, x, w, u) 9h(t, x, p, v) 

L ity = > Yi(°> x ) = a i( x ). xCG, P t y = — 


dz[ 


9 z i 


xC r, 


(4.9) 


M it z = - 


3 H ( t, x, w, u) n d 9 H(t, x, w,u) 

+ E “zr — + Pi( t - x )> 


9Yi k = i dx k 9y 


lx k 


3 h(t, x, p, v) n 3 H(t, x, w, u) 

z i ( T, x) = - Qi(x), Qiz — + £ — X k (x)l 


(4.10) 


9 Yi 9y 


lx k 


- Yt (t. x ). x c r. 


62 


Using the formula (4. 2) one can easily establish for arbitrary twice 
sectionally, continuously differentiable functions y^(t, x) and z^(t, x) 
that the Ostrogradskiy -Green formula 


m i r 

E JJ ( z x-itv + 

i = i o G 


y^M^z) dx dt 



-yj Qiz) d<r dt 


f Vi z i 

T 

dx 

G 

t = 0 


is valid. Let oo(r, x) = (u(t, x), v(t, x)) be a certain permissible control 
and y(t, x) and z(t,x) be the corresponding solution of the boundary 
problems (4. 9) and (4. 10). We will say that the permissible control 
Go(t, x) satisfies the maximum condition if 

H (t, x, w(t, x), u(t, x)) ((=)) sup H (t, x, w(t, x), u), xCG, O^t^T, 

uCU 


h(t, x, p(t, x), v(t, x)) (=) sup h(t, x, p(t, x), v), x£r, O^t^T. 

vCV 

where the symbol ((=)) indicates an equality which is valid everywhere 
within the domain C (0^ t ^ T, xCG) with the exception that there may 
exist points located on a finite number of n -dimensional surfaces, the 
( n -f 1 ) -dimensional volume of which is equal to zero. The symbol (=) 
is defined in an analogous way except that instead of n and G one 
should take n-1 and IT, respectively. The condition of minimum is 
defined in an analogous manner. 

THEOREM 7 (the maximum principle). For the permissible control 
cj(t, x) =(u(t, x), v(t, x)) to be min-op timal ( max -optimal) according to 
S, it is necessary that it satisfies the condition of maximum (minimum). 

This theorem, although it coes not supply the sufficient conditions 
of optimality, may serve as a practical instrument for the calculation 
of the optimum controls and the corresponding solution of the boundary 
problem (4. 4)-(4. 5). One can become convinced in this fact by 
repeating the deliberations carried out in Section I. 


2. The Formula for the Increment of the Functional S. The Proof of Theorem 7 

Let (t, x) be an arbitrary permissible control and y(t, x) and 
z(t, x) the corresponding solution of the boundary problems (4. 9) and 
(4. 10). Then 
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I = J £ ZiL it y - H( t, x, w, u (t, x)) dx dt 
C [i = 1 


/ 


c r 


m 


£ z i P iV “ h ( t > x > P> v ( tj x )) 


1 = i 


dcr = 0 , 


where 


C = (0^T, xCG), <r = (0^t^T f xCF), 


We take a certain permissible increment Ago = (Au, Av) of the control 
co(t,x) and denote by y-f-Ay and z+Az the solutions of the same problems 
(4. 9) and (4. 10) but corresponding to the control go -f Ago. Then 


m 


AI = I j^w + Aw, GO + AgoJ - I( w, GO ) ~J £ (Az^L^Ay + Az^^y + zjL^Ay) 

C * i=i 

— H(t, x, w + Aw, u + Au) - H(t, x, w, u)j j dx dt + /if (Az^Ay (4.12) 

J <r | i =i 


+ AziPj_y + z^Ay) - 


h(t,x,w+Aw, v+Av) -h(t.x w. 


i 


do - ~ 0 , 


while functions Ay^ and z-, i = 1 , . . . , m form the solution for the 
respective boundary problems 


8H ( t, x, w, u) 

L Lt Ay=A — , Ay^O.x) =0, xCG, 


3z| 


9h(t, x, p, v) 

PjAy = A , xC T , 


(4.13) 


dz^ 


M it Az= -A 


9H(t, x, w, u) d / 9H(t,x, w, u) 

: + £ — [ A 

9Vi k = i ox k 


8y ix k 


Az i (T,x)=0, x(;G, Qj t Az = A 


9h ( t, x, p, v) 


9yi 


(4.14) 


n 9H ( t, x, w, u) 



where 


aH 9H(t, x, w + Aw, u + Au) 9H(t, x, w, u) 


3 w k 3 w k 

ah ah(t, x, p f Ap, V + Av ) 

9 Pk 9Pk 


9w k 

ah ( t, x, p, v) 
3Pk 


Equality (4. 12) is now transformed by means of the formula (4. 1 1). 
Since the functions Ay and Az are the solutions of the respective boundary 
problems (4. 13) and (4. 14), we have 


Y f Az^LitAydxdt +J* AziP^Aydcr 

i = i LC o’ 


m 


£ J 

i= i ( C 


BH ( t, x, w, u) 

A A Yi 

3yi 


n d / 3H(t, x, w, u) 

ic — ( a 1 ^yi 


k = l dx k \ 


3Yi 


dx cit 


/ 


Bh (t, x, p, v) 

3>4 


n 9H ( t, x, w, v) 

+£ A x k<*> 

\r — t 


rlV; 


Ayjd(r I 


E | / 

i = i (c 


aH ( t, x, w, u) 

A Ayi 

3Yi 


n 9H ( t, x, w, u) 

+ Y A A y 


k = i 


3y 


ix k 


ix k 


dx dt 


ah(t, x, p, v) ) 

+ [ A Ayi do- [ • 

<r 3yi ' 


On the other hand 
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Y J Az^L^Ay dx clt +J Az^PjAydcr = J A 

1 = 1 lc o- J 1 =i[c 8zi 


A z^dx dt 


3h 

+ I A Az; dcr 

dz{ 


Consequently , 


V* f AziLqtAy ox dt +f Az^P^Aydcr 
i = l LC a 


r N - 3H 

E / * — 


L i = 1 C 9w i 


Awj dx dt 


(4.1 


2m - 3h 

+ E J A 9^ A Pi d(r 

i = l or 


where N = 2m -f nm - the dimensionality of the vector w. 
In an analogous manner we find 


S [/ Az i L ity dxdt +/^z i Piyd ( r 

i = 1 1 C, a- 


m 


=E 

i = l 


c 8H 

I Az;dx dt 

c azi 


c an 
+ J a z i 

c r 


A z^dcr 


(4.1 


m 


E 

i= l 


J* zjLjj-Ay dx dt +J z^PjAy dcr 
C or 


m 


E 

i = l 


/ Q i(x) Ayi(T, x) dx 
G 


+ / Pi(t, x) A yi (t, x) dxdt +J* Yi(t, x)Ayi(t, x) dcr 



The first sum on the right-hand side of the equality (4. 17) represents the 
increment AS of the functional (4. 6) during the transition from the 
control w(t, x) to the control <a(t, x) +Aco. Consequently, from the 
relationships (4. 12), (4. 15), (4. 16), and (4. 17) it follows that 


AS = 


- f H(t, x, \v 1 Aw, u + Au) - H(t, x, w, u) 

C L 

N / 9H 1 9H \ ] J 

-V ( — — + — A — 1 Aw; dxdt - / h(t, x, 

9w i 2 dwi / J J [ 

z ™ Id h 1 3h \ A ] j 

, x, p, v) (— — +— A— Api dcr . 

i=i \ 9 Pi 2 9 P X / J 


p + Ap, v 4 Av) 


- h (t 


Applying to the functions h, H, 9H/3w- and dh/dpi the Taylor 
formula and keeping within the expansion only the terms of the second 
order, one obtains, as it was done in £ 1, 


AS 


-J f H(t, x, w, u + Au) - H(t, x, w, u) 
C 1 


dx dt 



h ( t, x, p, v 4 Av) 


h(t, x, p, v) do- - r) , 


where n = ni + 


(4.18) 


hi 


J_ y. C /3H(t, X, w, u 

2 i = i c \ 9w i 


4 Au) 9H(t, X, w, u) 


Aw| dx dt 


+ 7 


1 f /dh(t, x, p, v 4 Av) 9h(t, x, p, v)' 


(4.19) 


2m - / 

Ef 

1=10- \ 


9 Pi 


9pi 


A pi do- , 
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l( f [9 2 H( t, X, w + 0iA w, u + Au) 

n2= 2. J — 


i,k=i C 


0 Wj^Wk 


0 2 H(t, x, w + 0 2 Aw, v + Av) 
3wj9w^ 


Aw^Aw^ dx dt 


2 m r 

£ / 


i,k=i 


9 2 h(t, x, p + 0 3 Ap, v + A v) 
0pi0pk 


9 2 h(t, x, p + 0 4 Ap, v + Av) 

0Pi0pk 


A Pi A Pk d(r 


To get the necessary estimates of the residual term r\ in formula 
(4. 18), we reduce the boundary problem to a system of integro- 
differential equations (see, for instance, work of Zagorskiy 34 , pp. 
90-96): 


r t dri 

Ay(t, x)= J / K„ (t,x f r, £)A — d£ 

0 G 


// K 12 (t, X, T.O tp (t, £) d^o- dT, x£G, 


an r r an 

*(t,X)= - A— + J J K 21 (t, X,t,£)A — d£ 
3 Z 0 G 02 


ff K 22 (t, X,r,e)4i (r.Odgo- dr , XCT, 

0 r 


where 


9H / 


9H \ 9h / 9 h 9h 

9z m / 9z \ 9Zi 9z m 


K ik - the Green type matrix. Inserting the value ip X) defined in 
(4. Zl) in the right-hand side of the same relationship and repeating 
successively a few times the same operation we obtain 
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(4.22) 


3h * f 8H 

4> (t, X) = -A + J f K n ( t, X, t ,r) ) A dt] dT 

9 Z o 6 •• 9z 


t 

+// K n (t, X, t,t|) iJj ( T,r|) d^cr dT 
o l' 

i r n_l i ah 

~//E K(t.X,T,T 1 )A-d T) adT I 

0 r i“0 


where 


t 

K n (t,x, T,T,)=K n _ 1 (t l X.T.T,) P) K 0 (a, P,T,r|) dptr da, 

r r 

t 

K n (t, X,t,t|) = // K n 1 (t, X,a, P) K° (a, p,T,ri) dgo- da, 
t r 

Ko=K 21 , K°=:K Z£ , n = 1 , 2, 

The number n is chosen so large that the kernel K n becomes 
bounded. This may be done because of the known estimates of the 
Green matrix and of its derivatives (see, for instance, work of 
Zagorskiy 34 , p. 92). Then from (4. 22) we get 


w( t) M Pj w ( t) dT ♦// Q n ( t > T > r i) 
0 o G 


3H 

9z 


d r\ dT 


t 

+/ / R n( t > T ’ r l) 


0 I 



ah 

+ r 

Ia — 

r 

| az 



d^cr dT 


(4.23) 


where P - a definite positive constant, 


(t) = X) | cl cr , Q »=/ 

r r 


R n : / £ | K M l » X.t.t,) 


d x (r . 
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We introduce the notation 


L ^ 

w k(t)=/ w k _ 1 (T)dT, w 0 (t)=w(t), Q n k = / Q nk-1 (t ' T > 1 l) dT > 

0 T 

t 

Qno = Qn> ■^•nk ~f ^nk -i (t.T.ri) dT, R m = R n (t,T,r|)+ 1, (4.24) 


V - 9 ^ 

XV — u J j | • • 


Integrating the inequality (4. 23) successively, we find: 


U L 

wk(t)^PJw k (T) dT +J J Qnk( t . T .'n) 


o G 
t 


♦// Rnk ( T,r| ) 

o r 


3H 

3z 

9h 


3z 


dr) dT 


do* dT . 


(4.25) 


Let us choose the number k so large that the functions Q n ^ and R n k 
be bounded for 0^ t^ t^ T, x^G and put 


Q(*> =// 


o G 
t 


3H 

9z 


f 

®l 

<i 

J 

r 

3 z 


max Q nk (9,T,r|) dr| dT , 
OS0-4 


max Rnk(0, T >ri) dcr. 
0- 0<t 


Then from the inequality (4. 25) for O<0s?t we get 

0 

w k (0)«ipj* w k (0)d0 +Q(t) + R(t). 
o 

We obtain from this using the known lemma (see work of Nemytskiy and 
Stepanon 30 , p. 19 ) that 


w k (0 )<- a( Q ( t) + R ( t)| 

for 0^0 - t^T and, consequently, 
w k (t) <- A| Q(t) +R(t)|, 
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where A is a definite positive constant. Thus it follows from the 
relationships (4. 23), (4. 24), and (4. 25) that 


w k< 


/ M 1 (t,T,Tl) 

aH 

A 

dr] +J Nx ( t, t,t] ) 

3h 

A 

dcr 

Lg 

3z 

r 

3z 

- 


, (4.26) 


where NLi and Ni are functions of the same type as Q n and R n , 
r e spectively. 

Since the number n was chosen sufficiently large, we find from 
the relationships (4. 22) and (4. 26) that 


h (t,x)p 


t r 

/ dT • 

f M 2 (t, T, X,T|) 

aH 

A 

0 k 

Jf 

3z 


drj 


+ / N Z (t, T.X.Tl) 

r 


ah 

A 

dz 


d T]°" 


(4.27) 


where M 2 and N 2 - the scalar functions of the Green type functions. 

The functions 3H/9z and dh/dz are continuous over t and twice 
continuously differentiable over other arguments. Consequently, for 
each permissible control (4. 20) may be differentiated also with respect 
to Xi, . . . , x n and by means of the above described method one can 
obtain the inequalities 


| A gi(t, x )j^/ M 3 (t, X,T, t)) £ Aus (T, r|)| dT dr\ 
C S= 1 

+ / N 3 ( tj X ’ T,T l) ^ | Av j ( Tj ^^1 d T,r ) 0 ' ' 

c r J = 1 


(4.28) 


where 


g = 



Ym. 


9Yi 

8 xj 


x£G, 0«t<T. 

9x n / 


Since the functions Az^(t, x), i = 1 , m, form the solution of 

the boundary problem (4. 14), we find in an analogous manner that 
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(4.29) 


) Az i (t, x)|^ J M 4 (t, x, t, Ti) £ |Au s (t,t|)| d-n dT 
C s = 1 


+ f N 4 (t, x,T, T|) J 2 (Avj {r,T|)| d TjT 1 o- , 

(T j “ 1 


x£ G, 0 t « T, i = 1 , . . . , m. 

Because of these inequalities, we find from (4. 19): 

M« b i/ £ | Au j(t»x) il/ M U (t, X, T, T] ) 2 

C j — i (c j=i 


Auj (t, r)) 




q 

Nn (t, X,T, T]) £ 

k= i 


^ V k d T>7]W ( dx dt 


q ( r 

+ B z f ^2 l Av i( t » x )| !/ Mn(t,x,T fTI )i: |Auj|dr|dT 

cr i =i ( C j =i 

+f Ni i ( t, x,T,T,) | A v p| d T)r j(r | dj-jjjcr , 

cr k = l ) 

where - positive constants, 

Mn = M 3 + M 4 , Nn = N 3 + N 4 . 

Since according to the condition 


3 2 H 

and 

3wj3 wj 


are bounded, then 


3 2 h 

3Pi9pj 


dr] dT 


(4.30) 
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m 


B 


/ / M n (t, x.T.q) E 
C C k =1 


Au k 


dT dr) 


tJ*N n (t, x,T,q) £ 


J = i 


AV J 


d T, n 0 " 


dx dt 


(4.31) 


4 B 4 r / M u (t, x,T,q) *E |AuJ dq dT 
a C k=i 1 1 


I / N n (t, X,T, T) ) E 

or j " 1 


Av 


d T,r| (r 


dt, xo- 


If one takes into account that the constant B may be chosen such that 


r 2 q 

J Au k (t,x) dxdtxB 2 , E J|Av:(t, 
k 1 C 1 1 j = 1 0- 1 


X) 


2 do- ' B 2 , 


then from the inequalities (4. 30) and (4. 31) it follows that for the 
residual term q in the formula (4. 18) one has an estimate: 


T / JpU.x.T.nJE |Au k 

1 C ( C k = i 1 


(It dr) 


r J Q ( f. x » T > [ |) E 

O' j =1 


r 


U T,q ff j C1Xnt 


/!/ P(t, x, t, q) E I A u k I (iTdq 


(4.32) 


(T ( C 


k =: 1 


/ Q ( t, X, T, q ) £ I A Vj 


J 


d T,q°( d t,x^ - 


where the functions P and Q are of the same type as Mu and . 


Formula (4. 18) and inequality (4.32) are analogous to the 
corresponding relationships in Section I . Consequently, the proof 
of Theorem 7 conincides almost literally to the proof of Theorem 1. 


If the system of equations (4.4) is linear and has the form 


m, 


(4.33) 


m 

L it y = X) d ik (t> x) Yk + f i ( v ), i = 

k = 1 


while the boundary conditions may be represented in the form 


m 

Pi(t,x)y=£ c ik (t, x)y k 4 cp.(v), 

k = l 

(4.34) 

X c r, y ( 0, x) = a. (x), xCG, 


then the following theorem is true: 

THEOREM 8. If to each permissible control corresponds a unique 
solution of the boundary problem (4. 33)-(4. 34), then, for the control 
oo(t, x) = (u(t, x), v(t,::)) to be min -.optimal (max-optimal) according to 
the functional (4. 6), it is necessary and sufficient that it satisfies the 
conditions of maximum (minimum). 

The proof of this theorem follows directly from the fact that in 
the case under consideration the formula (4. 18) for the increment of 
the functional S takes the form 


AS = - J* | H ( t, x, w, u + Au) — H(t,x, w, u)| dx dt 
C 

(4.35) 

- J | h ( t, x, p, v -f Av) - h (t, x, p, v) | dcr . 
cr 

3. Problems with other Optimality Criteria 

The just obtained result may be applied to the solution of 
problems of optimum control with other optimality criteria. 

Let, for instance, the control process be described by the boundary 
problem (4. 4) -(4. 5) for which the domain G is a rectangle O^x^X^ 
and let one choose as the optimality criterion the functional 

T XiX 2 

* Iff fo (t, x, y, y X) u) dx 2 dxidt. (4.36) 

0 0 0 

We introduce an auxiliary variable y 0 by means of the relationship 
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= fo (t, X, y, y x> u), y 0 (xi, x 2 , 0) = y 0 (x!, 0, t) = y 0 (0,x 2 , t)= 0. 

9x! 9x 2 9 1 


Then the problem reduces to the evaluation of the minimum of the 
functional S = yo ( , X 2 , T). We form the function H: 

m 

H ( t, x, w, u) =2^ z i f i (t, x, y, y x , u) + z 0 f 0 (t, x, y, y x , u). 
i = l 


The function z^(t, x) is defined by means of the equations 

9 3 z 0 

9x 1 9x 2 at 

and additional conditions (see formulas (2. 19) and (4. 10)): 


9H 


M it z = 


9H 


9 yi k = l dx k Uy 


i = 1. 


ixk, 


= 0 


9h(t, x, p, v) 3 3H(t, x, w, u) 

Q it z = + E X k( x )> x cr , 

9yi k = i 9yix k 


0 Z 0 0 Z Q 

= 0 for xi = Xi , x 2 = X 2 , = 0 for t = T, x 2 = X 2 ; 

9t 0xi 


0Zq 0 2 Z O 0 2 Z O 

-= 0 for t = T, Xl = Xi , — _ 0 for t = T, — = 0 

0 X 2 0 Xj[ 9 X 2 0 Xi 0 t 


d z z 0 

for x 2 = X 2 , = 0 for x x = Xi , z 0 (X x , X 2 , T) = -1, 

0 X 2 0 1 


z i (x lf x 2 , T) = 0, i = 1, . . • , m. 

In this manner, z 0 (xi , x 2 , t) = - 1 , and the function H take s form 
H = H (t, x, w, u) - f 0 ( t, x, y, y x , u) , 

where 

m 

H = 23 z ° f ° ( ' t ’ x > y> yx. u ) ■ 

i = i 


Let us study the functional 
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I = 


T XjX 2 

Iff 

0 0 0 


£ 3 3z o 

Lj Z i-L'i t y + z o T — I — IT ~ H ( t, x, w, u) 


1 = i 


9xi8x 2 8t 


dx 2 dx!clt 


♦// 

o r 


^ z i p iY ~ h (t. P. v ) 


1 = 1 


clcrdt = Ii + I 2 , 


where 


TX,X Z 


■. -/// 


0 0 0 


m 


X z i L ity -H(t, x, w, u) 


1 = 1 


dx z dxi dt 


♦// 

o r 


£ ZiPi(t, x) y -h(t, x, p, v) 
i = 1 


dcrdt , 


TX,X Z 

12 = / / / 

0 0 0 


9 3 yo 


8xi 9 x z 8x 3 


- f 0 (t, x, y, y x , u) 


dx z dxidt 


By transforming the integrals I L and I z in the same manner as done in 
Sections I and IV, we obtain a formula for the increment of the 
functional (4. 36) in the following form: 

T XiX 2 

AS = -/// j H ( t, x, w, u + Au) — H(t, x, w, u)| dx z dxi dt 
ooo 

T 

-J J | h(t, x, p, v + Av) - h(t, x, p, v)| da dt - r\ , 
o r 

where the residual term r\ is defined by formulas analogous to (4. 19 ). 

Consequently, the necessary conditions of optimality for the 
problem under investigation may be formulated in the form of Theorem 
7 where in the conditions for the maximum (minimum) of the function 
H is substituted by H. 

Using the results of Section I, one can analogously investigate 
other problems of optimum control processes where for the optimality 
criteria one utilizes various nonlinear functions. In particular, the 
results obtained may be applied to the study of problems investigated 
by Bellman and Osborn 5 and by Bulkovskiy Lerner 6 . 
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4. Optimum Problems in the Theory of Elliptical Systems 


The control problems are analogous to those which were 
investigated above and are encountered during the study of diffusion 
processes 3 ’ 8 . Nevertheless, one must investigate boundary problems 
for elliptical equations. Problems of such a type are encountered 
during the study of optimum thermal and electrical fields in various 
power devices. 

At this point we will briefly outline the formulation of the 
maximum problem for electrical systems and derive the formula for 
the increment of the functional by means of which one finds the 
optimality conditions. 

Thus, let us deal with an elliptical system of equations 

L,y = f(x, y, y x , u), x = (xj, . . . , x n )£ G , (4.37) 

where the operator L is defined by the formula (4. 1) and the function 
f = (fi , . . . , f m ) is twice continuously differentiable over the totality 
of all its arguments. The control parameter u takes the values from 
a bounded domain U (closed or open) of an r -dimensional euclidian 
space. 


Let further the function y(x) satisfy the boundary conditions 

Pi(x) y = Ti(x, y, v), i = 1, . . . , m, xCT , (4.38) 

where CHq satisfies the same conditions as fp and the parameter v takes 
values of a bounded domain V of the q -dimens ional euclidian space. 

The permissible control oo(x) - (u(x), v(x)) is defined in the same 
manner as in Part 1, and we assume that over the discontinuity surfaces 
of the control the desired function satisfies certain smoothness con- 
ditions 35 . We assume that we impose onto the known functions of the 
boundary problem, in addition to the above mentioned conditions, some 
additional limitations under which for each permissible control there 
exists a unique solution of the particular problem. 

We formulate the problem: among all the permissible controls 

we determine the control co(x) (if it exists) such that the corresponding 
solution y(x) of the boundary problem (4. 37) -(4. 38) realizes the 
minimum (maximum) of the functional 
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(4.39) 


m 

s=E 

i = 1 


J a i (x)y i (x)dx + f Yi(x)Yi( x ) d o- 

LG r 


where ctj,(x) and Yj.( x ) are gi- ven continuous functions. 

We introduce the functions H z^f^ and h z^cpp The 

function z^(x) is defined as the solution of the boundary problem 


M { z 


£H 

9yi hi 


d / du 
dxkl 9yix k 


Qi(x), x£G , 


(4.40) 


9h n aH 

Qi z = T) X k( x ) + \i( x ) > x cr (4.41) 

9Vi k=i 9y ixk 


(the definition of the operators Mj and can be found at the beginning 
of the paragraph). If it appears that the right-hand side of Equation 
(4. 40) contains the derivatives of v Xi (x), . . . , v Xn (x), then one demands 
from the permissible controls that they have sectionally continuous 
derivatives with sufficiently smooth discontinuity boundaries. Then 
the boundary problem (4. 40)-(4. 41) for each permissible control has 
a unique solution. 

In the same manner as it was applied above, one can obtain a 
formula for the increment of the functional (4. 39) in the following form: 


AS = -/ j H(x, w, u + Au) - H ( x, w, u)|dx 
G 

-J| h(x, p, v + Av) - h (x, p, v)j dcr - q 

T 


(4.42) 


where the residual term q is determined using formulas analogous to 
(4. 19). If the boundary problem (4. 37 ) -( 4. 38) is linear, then rj = 0, and 
consequently, we have as valid: 


THEOREM 9. For the permissible control to be locally min- 
optimal ( max -optimal) according to the functional (4. 39) in the linear 
boundary problem (4. 37)-(4. 38) (functions f^ and ^ are linear in y 
and y x ) it is necessary and sufficient that this control satisfies the 
conditions of maximum (minimum). 
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In conclusion we note that the analogous problem (with analogous 
results) may be investigation also for the system of hyperbolic equations 
with initial and boundary conditions. 
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Section V. SOME PROBLEMS OF INVARIANCE THEORY 


Let the control process be described by a system of equations with 
partial derivatives 

Az = f (xi x^, z, u), (5,1) 

where A is a linear differential operator of parabolic, elliptic, or 
hyperbolic type, z = (zi , . . , , z m ) is a vector characterizing the state 
of the system under control, and u is a vector characterizing the exter- 
nal interaction. Let also be given additional conditions which may 
contain vector v defining the external interaction on the system. We 
assume that the vector co = (u, v) is subjected to the same conditions as 
the permissible control in the problems of optimum control discussed 
above while the additional conditions are such that to each vector co 
corresponds a unique solution of Equations (5. 1) with the same additional 
cond itions . 

Let, in addition, be given a certain functional I [?] defined over the 
solutions of Equations (5. 1). The basic problem of the invariance theory 
is to find conditions for which the functional I does not depend on the 

1 7 

external interactions. In the paper by Rozonoer it was shown that 
the invariance problem may be studied by the methods of the variational 
calculus in the case when the control processes are described by ordi- 
nary differential equations. Analogously, one may investigate the 
invariance problem also for the system with distributed parameters. 

We investigated the control system whose behavior is described by 
the boundary problem (4.9) with certain smoothness conditions imposed 
on the discontinuity surfaces of the function u(t, x). We assume that 
with each permissible vector go (t, x) - (u)t, x), v(t, x)) is associated 
a respective unique solution of the same boundary problem and 

m 

fi (t, X, y, y x , u) d ik( t ’ x ^k + 8i< t > x ) u > 

k = l 

(5.2) 

CP i ( t, x, y, v) - p i (t, x)v, 

where, for the purpose of simplification of subsequent formulas, u and 
v are viewed as scalar quantities. 

As the functional I we utilized the expression (4. 6) in which the 
time T and the domain G are viewed as fixed. The "adjoint" boundary 
problem (4. 10) has in this case the form 


80 


m 

M u z r -23 r 1 ki z k 1 V 1 ’ x) ’ z i< T » x) x e g, 

k l 

(5. 3) 

Qjz - Yi(t, x), xc r. 


The formula (4. 3 5) for the increment of the* functional (4.6) takes the 
for rn 


AS 


/ m 

-/*«(£ Oj/jj dxdt -/* Av 23 Pi z i 

C 1 1 cr 


d o' 


Consequently, if 


m 


23fii( l > x ') 7 i( f > x ) = 0, > C G 


1 -- 1 


X) Pi(G x ) z l('> x ) = °> xCI’, O' t< T, 


L - 1 


(5.4) 


the functional S does not depend on the external perturbation co(t, x). 

By using the method of proving thc^ opposite, it is easy to establish 
that these conditions arc' also necessary for the* functional not to depend 
on u> (see, for instance, work of Rozonoer lf ). For the verification of 
the c ondition ( 5 . 4 ) one must find a solution of the boundary problem (5, 3). 
However, for the special case* presented below, one can establish the 
necessary and sufficient invariance* conditions expressed through the 
coefficients of the* equations of the* boundary problem (4. 9). 


Let the control process be* described by the equations 


m / 9u. n d y. 

Ly i E <' ik v k * RjU pi - D Ojk ^ 

k -- i ' j , k l J 

d ik> Si - tonst ) 
with the additional conditions 

Yj (0, x) = a i (x), x C G, Pyj -■ p (t, x), xCT; 


3 Yi_ 
9x k 


(5. 5) 


(5. 6) 
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here P is a linear differential operator defined over the boundary T, 
where the differentiation is carried out in a direction which is towards 
the outside relative to G. 


Let us study the functional 


S = 


m r T 

T J a i( x )Yi(T, x) dx + if 'Y i (t, xjy^t, x) dcrdt 

i= i Lg or 


(5. 7) 


Then the functions z^(t, x), entering (5.4), are determined from the 
equations 



with the additional conditions 

z i(T , x) - - a i(x), x C G; 

Q(t, x) z i = - Yi(t, x), x C r, O v t T, 


(5.8) 


(5. 9) 


where the operator Q is defined according to (4. 3) as being conjugate 
to P. 


We denote by D and D* the matrix of the coefficients d.^. and their 
conjugate matrix, respectively, and by (r, s) - the scalar product of 
the vectors r and s. Since in the case under investigation the boundary 
conditions (5.6) do not contain v, the conditions (5.4) may be written 
in the form 


R(t, x) - (z, g) - 0, xC G, 0 < t < T. 

Applying to this equality the' operator M and taking into account 
that z satisfies the system of equations (5. 8), we get 

MR(t, x) - (Mz, g) - - (D-z, g) - (z, Dg) =. 0, 


Analogously, we find that 



M k R(t, x) = ( - 1)^ (z, D^g) = 0, k-0, 1 ,..., m - 1 . (5. 10) 

From this, because of the conditions (5. 9), it follows that 

M k R(T, x) = -(-l) k (a (x), D k g) = 0, k = 0, 1 , . . . , m - 1. (5.11) 

Putting in Equations (5. 10) x C rand applying the operator Q, we obtain 
taking into account the conditions (5. 9): 

QM k R(t, x)= (-l) k (Qz, D k g) = -(-l) k ( Y (t, x), D k g)= 0, 

(5. 12) 

k = 0, 1 ,..., m - 1 . 


The conditions (5. 11) and (5. 1Z) are necessary for the invariance 
of the functional (5. 7) relative to the external interaction u in the boun- 
dary problem (5, 5) - (5. 6). Let us show that these conditions are also 
sufficient. 


Since according to the condition even one of the vectors a = 

(ai, . . . , a m ) and y - (yx, . . . , y m ) differs from the zero vector, there 
exist numbers \ Q , . . . , \ m _ \ such that 


m-i 

Z ^r d g = °- 

k = o 

Multiplying the k-th equation (5. 10) by (-1) \ k and summing over all 
k, we obtain: 

m-i 

£ (- 1) \ k M R(t, x) = 0, xCG, 0 t • T. 

k- o 


Introducing the notations 

R k - M R, k = 0, 1 , . . . , m-2, 

we obtain from this equation and the conditions (5. 10) and (5. 11) a homo- 
geneous boundary problem for the determination of R^: 
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(- 1 ) 


m- i 


, m- 2 


m 


_ i MR m _ 2 + ( - 1 ) ^•m - 2 R m - 2 + . . . - ^iRi + X- o R o - 0 


R m _ 2 (T,x) = 0, x£G; QR m _ 2 ( t, x) = 0, x C r, 


MR m-3 ' R m- 2 = °> R m-3^ T > x ) = °> x £ G; QR m _ 3 (t, x) = 0, xC T, 


MR 0 - Ri = 0, R 0 (T , x) = 0, x CG; QR 0 (t,x) = 0, x £ r. (5. 13) 

Since we assume that the coefficients of the operators M and Q are 

sufficiently smooth, the boundary problem (5. 13) has only a trivial 

3 4 

solution (see, for instance, work of Zagorskey , pp. 97- 103): 

R m _ 2 (t, x ) = Rm-s^’ x ) = • • • = R o (t> x )- °> 
and from this follows the validity of the' condition 
m 

Z £i(t, x ) z i(t. x ) = 0. x£G. 
i - 1 

This proves: 

THEOREM 10. Eor the invariance of the functional (5.7) relative 
to the external interaction in the boundary problem (5. 5)- (5. 6) to be 
true, it is necessary and sufficient that the conditions 

(a (x), D k g) =0, x £ G; 

(Y (t, x ), D k g) = 0, x £ r; 

0 : t^ T , k = 0 , 1, m - 1. 

be fulfilled. 


From the method for the proof of this theorem is clear that the 
analogous results may be obtained for boundary problems which were 
investigated in Section I. In particular, for the boundary problem (1. 33) 
one has to utilize the formula (1. 35) for the increment of the functional 
(1.3) with the function u^ defined by means of the boundary problem 
(1.34). 
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